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ABSTRACT 

We provide evidence for a correlation between the presence of giant clumps and the occurrence of ac- 
tive galactic nuclei (AGN) in disk galaxies. Giant clumps of 10 8-9 M© arise from violent gravitational 
instability in gas-rich galaxies, and it has been proposed that this instability could feed supermassive 
black holes (BH). We use emission line diagnostics to compare a sample of 14 clumpy (unstable) 
disks and a sample of 13 smoother (stable) disks at redshift z~0.7. The majority of clumpy disks 
in our sample have a high probability of containing AGN. Their [Om]A5007 emission line is strongly 
excited, inconsistent with low-metallicity star formation alone. [Ncin]A3869 excitation is also higher. 
Stable disks rarely have such properties. Stacking ultra sensitive Chandra observations (4 Ms) reveals 
an X-ray excess in clumpy galaxies, which confirms the presence of AGN. The clumpy galaxies in 
our intcrmcdiatc-redshift sample have properties typical of gas-rich disk galaxies rather than mergers, 
being in particular on the Main Sequence of star formation. This suggests that our findings apply 
to the physically-similar and numerous gas-rich unstable disks at z > 1 . Using the observed [Olll] 
and X-ray luminosities, we conservatively estimate that AGN hosted by clumpy disks have typical 
bolometric luminosities of the order of a few 10 43 ergs -1 , BH growth rates tobh ~ 10 -2 M yr - , and 
that these AGN are substantially obscured in X-rays. This moderate- luminosity mode could provide 
a large fraction of today's BH mass with a high duty cycle (>10%), accretion bursts with higher 
luminosities being possible over shorter phases. Violent instabilities at high redshift (giant clumps) 
are a much more efficient driver of BH growth than the weak instabilities in nearby spirals (bars), and 
the evolution of disk instabilities with mass and redshift could explain the simultaneous downsizing 
of star formation and of BH growth. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: ISM — galaxies: high-redshift - 
galaxies: formation 



1. INTRODUCTION 

The formation of Supermassive Black Holes (BH) 
at the center of galaxies appears to be connected to 
their stellar mass assembly history. Not only is to- 
day's BH mass close l y related to the stellar mass 
(jMagorrian et all 119981 : iMerritt fc Ferraresd l200l . but 
also the BH accretion history appe ars to follow the cos - 
mic star formation (SF) history (jBarger et al.l |2001[ ). 
The bulk of both SF activity and BH growth has 
long been attributed to mergers, especially at high 



redsh ift (e.g.. lElbaz fc Cesarskvl 120031 : iDi Matteo et al.1 
[2001 iHookins et al.ll2006l )" but there is increasing recog- 
nition that the SF history is not dominated by merger- 
induced starbursts, but by continuous SF in gas-rich star 
forming disks (e.g.. lElmegreen et all 120071 : iDaddi et al.l 
[2?MlRodighiero et al.ll201lTr Models indicate that mas- 
sive galaxies at high redshift obtain a large fraction of 
their gas from rapid cold accretion, most of which comes 
as smooth gas streams and very small mergers, not ma- 
jor mergers: the high rate of continuous gas infall main- 
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tains high fractio ns of cold gas in the disk and high star 
formation rates (IKeres et al.1 120061 [Brooks et aLl 120091: 



IDekel et al.|[2009l : IBouche et al.ll201ofr 

Hence, galax y mergers likely pr oduce the strongest 
starbursts (e.g.. lTacconi et aT1l2008f ) but they appear not 
to dominate the cosmic star formation history. Similarly, 
the most rapidly accr eting BHs in bright QSOs are prob- 
ably fed by mergers (|Riechers et al.l 2008; Walter et al.1 
l200l . but mergers do not necessarily dominate the 
global budget of BH growth : moderate AGNs with sub- 
quasar luminosities at z sa 1 — 3 mostly reside in nor- 
mally st ar-forming galaxies rath er than in starbursting 
mergers (jMullanev et al.ll2012al) . While it is observed, 
at least locally, that mergers and interac tions do trig- 
ger A GN compared to isolated galaxies (|Ellison et al.1 
1201 lh . mergers may not be the dominant feeding pro- 
cess, due for instance to the scarcity of major merg- 
ers compared to the general galaxy population. In- 
deed, X-ray AGN host galaxies do not show morpho- 
logical merger signatures more frequently than inactive 
galaxies, especially at high redshift and at least for the 
most frequent moderate-luminosity AGNs dGrogin et al.1 
120051 : iGabor et al.l 120091 : iCisternas et al.l 120111 ). The X- 
ray signal of AGNs coul d be attenuated by obscura- 
tion in gas-rich mergers (jHopkins et al.ll2006D . but the 
most recent studies with increasing X-ray sensitivity 
confirm the lack of any observable AGN-merger con- 
nection dSchawinski et al.1 120111 iMullanev et al.1 12012a! 
iKocevski et al.1 1201 If ). Non-merging galaxies also ap- 
pear to dominate the X-ray luminosity and BH accre- 
tion density, especially for moderate-luminosity AGN 
(|Georgakakis et ai]|2009D . All these observations suggest 
that high-redshift BH growth largely proceeds indepen- 
dently from major galaxy mergers. 



Could internal processes feed BHs and AGNs effi- 
ciently in high-redshift disk galaxies? Internal, "secu- 
lar" AGN feeding in nearby spiral galaxi es mostly involve 
bars (and embedded nuclear bars, e.g.. lEmsellem et al.1 
I200H ): this is a slow process, exp ected to provide only 
a sm all fraction of the BH mass (jHopkins fc Hernquistl 
120061 ). Yet, disk galaxies at high redshift are very differ- 
ent from nearby spirals. They are very gas-rich, strongly 
unstable, often dominated by giant clumps of gas and 
young stars - while local disk galaxies are dominated 
by bars and spiral arms, with only low-mass star form- 
ing complexes inside the spiral arms. The ubiquity 
of clump-dominated morphologies in high-redshift star- 
forming ga laxies has be e n poin t ed out at various wave- 
lengt h s bv iCowie et al.l (119961 ) : lElmegreen et~aTI (12004 
l2007l):lForster Schreiber et al.1 (j2006l 120091 ) : lGenzel et all 
(|2008f ): IGuo et al.l (120111 ). These giant clumps are 
not just low-mass patches of star formation: they 
have individual masses up to a few 10 8 ~ 9 each 
(lElmegreen fc Elmegreenl 120051: lElmegreen et~"aH l2009at 
iForster Schreiber et al.1 120061 I2011D . Most clumps are 
not incoming satellites, but form in-situ by grav- 
itational instabiliti es, as indicated by their p hoto- 
metric prop erties (lElmegreen et all 120071 l2009aD, gas 
kinematics dShapiro et al.1 120081: iBournaud et al.1 [200l 
Ivan Starkenburg et al.l 20081) . and est imates of Toomre 
parameters Q < 1 ( Genzel et al.ll201ll) . 

Giant clumps are actually the most striking feature 



tracing the global, violcnlQ gravitational instability in 
high-redshift disks. This instability results from very 
high gas fractions of th e or der of 50% as obser ved by 
IDaddi et al.1 (120081120101 ) and lTacconi et al.l (j2010f ). The- 
oretically, the high gas fractions and instability are ex- 
plaine d as the natural ou tcome of the rapid cosmic infall 
of ga s (IDekel et al.ll2009l (her eafter DSC09); lAgertz et al.1 
I2009t iCeverino et al.1 I2010D . while merger-dominated 
phase s produce stable systems (jBournaud fc Elmegreenl 
120091 ). The disks become more stable and clumpy struc- 
tures are replaced by spiral arms an d bars at lower red- 
shift, as the gas f raction decreases (jMartig et al.l 120121 : 
iKraliic et~a l. 2012) and the turbulent velocity dispersion 
of the gas also decreases (Kassin et al., in preparation). 

It has been recently proposed that disk instability in 
such high redshift conditions could efficiently feed BHs. 
The gaseous and stellar densities in the giant clumps 
are so high that black holes may start forming by run- 
away stellar collis ions and become the see ds for future su- 
permassive BHs (jElmegreen et a l. 2008a). Furthermore, 
the observed instability and the associated turbulent en- 
ergy dissipation trigger a gas inflow throughout the disk 
towards t he nucleus, via intense gravitational torquing 

" ' 20091 Krumholz fc Burkertl 
2010I ). The central inflow 



dGammid 120011: IDekel et al 
I2010t lElmegreen fc Burkert 



rate could feed a BH at a realistic growth rate over a 
few 10 8-9 yr, along with bulge growth (Bournaud et al. 
2011, hereafter Bll). This gas inflow toward a central 
BH can be triggered even when giant clumps are present 
in the outer disk, regardless of the clumps being long- 
lived against stellar feedback and able to migrate radi- 
ally (jKrumholz fc Dekelll2010l Bll). The cold gas infall 
maintaining the ins tability may persist even with power- 
ful A GN feedback (jDi Matteo et al.l 120111 iDubois et all 
I20T2T ). 

This mechanism could hardly be directly tested in z « 2 
clumpy disk galaxies using current facilities, though. The 
associated AGNs should most of the time have modest 
intrinsic X-ray luminosities, and would be strongly ob- 
scured by the interstellar medium in these high-redshift 
disks, which are gas-rich and geometrically thick (Bll), 
in addition to potential obscuration by the AGN torus 
itself. Instead, narrow emission lines from gas surround- 
ing the AGN could be a better tracer of moderate BH 
activity in gas-rich galaxies: the associated emission re- 
gion is more extended spatially, and the existence of un- 
obscured sight lines is more likely - the use of narrow 
line diagnostics to probe mode r ate A GNs was illustrated, 
for instance, bv iWright et al.l (120101) . However, the ap- 
plication of these diagnostics to primordial galax ies at 
z ~ 2 remains uncertain (e.g.. iShaplev et aLll2005l) . and 
the main emission lines are redshiftcd in the near-infrared 
where available data arc limited to few systems with low 
signal-to-noise spectra. The SINS near-infrared survey 
gathered near-infrared spectra with high signal-to-noise 
ratios for a sample of z ~ 2 star-forming galaxies, but 
it is mostly targeted toward galaxies that are thought to 

We call "violent" the ring, clump, and axisymmetric insta- 
bilities in high-redshift disks, in contrast to the weaker, non- 
axisymmctric modes (bars and spiral arms) that dominate low- 
rcdshift disks. In the violent high-redshift instability, gravitational 
torques arc considerably stronger and the radial gas inflow is much 
faster (Bll). 
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be fr ee from AGN, on purpose (iFdrstcr Schrcibc r et al.l 
I2009D - although it does contain at l east one example o f 
AGN in a rotating clumpy disk (e.g. JGenzel et al.1l2006l ). 

To get around these problems, we study in this pa- 
per the presence of AGNs in clumpy disk galaxies in the 
intermediate redshift range 0.4 < z < 1.0. An advan- 
tage of this intermediate redshift range is that emission 
line diagnostics can be used with higher confidence than 
at z ~ 2, because the mass-metallicity relation evolves 
at a lower rate at z < 1 than at higher rcdshifts (e.g., 
iSavaglio et al.l[2~005| ). Another advantage is that in this 
redshift range, stable disks more similar to nearby spirals 
are also common, making it possible to directly compare 
violently unstable, clumpy disks to smoother, more sta- 
ble systems. Clumpy galaxies in this redshift range are 
fully representative for the global process of violent disk 
instability at high redshift (as further demonstrated with 
our own sample in Section 2). Hence, focusing on inter- 
mediate redshifts enables us to use reliable AGN diag- 
nostics, apply them to samples of clumpy unstable disks 
galaxies, and compare to more stable spiral-like galaxies 
at similar redshifts. 

In Section 2, we build a sample of Clumpy galaxies and 
a reference sample of Stable disks in the GOODS3-South 
field. The dumpiness is measured both visually and by 
computer; the samples are about mass- and redshift- 
matched. In Section 3, we use narrow emission line di- 
agnostics, in particul ar the Mass-Excitati on (MEx) diag- 
nostic introduced bv I Juneau et al.1 (|2011l . hereafter Jll), 
which is useable with optical spectra up to redshift one, 
statistically calibrated, and robustly tested against X- 
ray AGN selections up to z«l in Jll. In Section 4, we 
perform X-ray stacking using the deepest 4 Ms Chandra 
data. The inferred AGN luminosities and BH accretion 
rates are discussed in Section 5. 

Throughout the paper, we assume f2 m = 0.3, £Ia = 
0.7, H = 70 km s" 1 . [Om] denotes [Om]A5007, [Nem] 
denotes [Nein]A3869, and [On] denotes [On]AA3726,3729, 
the [On] flux being the total flux of the doublet. Line 
ratios such as [Om]/H/3 denotes a flux ratio, except for 
[Oll]/H/3 which is an equivalent width ratio. 

2. CLUMPY DISKS AND STABLE DISKS: SAMPLE 
SELECTION, CLASSIFICATION, AND ANALYSIS 

Violent instabilities and giant clumps are most fre- 
quent above redshift one. Here we select intermcdiatc- 
rcdshift clumpy systems at z ~ 0.7, in order to uti- 
lize robust diagnostics that can distinguish the signatures 
of AGNs from SF, and compare to a control sample of 
more stable disks that are absent from z > 1 datasets. 
In studies covering a large redshift range, the morpho- 
logical and photometric properties of clumps in galaxies 
at intermediate redshift are not found to be very differ- 
ent f r om tho se of z ~ 2 clumpy disks (|Elmegreen et all 
120071 l2009al) . The only difference is that violent insta- 
bilities persist at lower redshift only in moderate-mass 
galaxies. At redshift 2, clumpy disks are observed even 
in systems with stellar masses of several 10 11 M (7 1 (e.g., 
iGenzel et all [20061 ; lF5rster Schreiber et al.l I2009T) . Be- 
low redshift one, clumpy galaxies typically have baryonic 
masses of the order of a few 10 M Q (e.g., Elmegreen 
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et al. 2009a, Puech 2010). The fact that clumpy mor- 
phologies persist down to lower redshift for lower mass 
galaxies is consistent with the theoretical framework in 
which violent disk instabilities arise in gas-rich disks fed 
by cold streams. The persistence of the instability pri- 
marily re quires the preservation of a hig h gas fraction in 
the disk feournaud fc Elmegreen! [20091 . DSC09). Cold 
accretion shoul d persist down to lower r edshift for lower- 
mass systems (jDekel fc Birn boim 2006(). More massive 
gala xies form their star s and consume their gas earlier 
(e.g. Uuneau et al.l 2005 ). indeed their gas fractions are 



lower (jKannappar 



2004T ) and are unlikely to support vi- 



olent instabilities. More massive galaxies also build a 
stellar sphero id more rapidly, which also stabilizes the 
disk fDSC09. iMartig et alJl200iCacciato et al.ll2012ft . 

Hence, the fact that clumpy galaxies at intermediate 
redshift are mostly moderate-mass systems (a few 10 10 
Mq) while such morphologies are also found in more mas- 
sive galaxies at z ~ 2, does not indicate a different ori- 
gin or a different dynamical evolution. This is actually 
a natural property of the global process of violent disk 
instabilities in gas-rich galaxies fed by cold gas flows. As 
lower-redshift systems tend to have lower gas fractions 
and become more stable, an external tidal field in distant 
interactio ns may sometimes be needed to trig ger the in- 
stability (|Di Matteo et all 120081: iPuechl 120101 ) . but even 
so the dynamical process remains a violent instability in 
a gas-rich disk. We will show in the following that our 
sample of z ~ 0.7 clumpy galaxies are gas-rich systems 
on the Main Sequence of star formation, as expected for 
clumpy disks at higher redshift. 

2.1. Strategy and Datasets 

We aimed at obtaining a sample of Clumpy unstable 
galaxies at z ~ 0.4—1, and a comparison sample of Stable 
disk galaxies with only weak non-axisymmetric instabil- 
ities (bars and spiral arms). Having two such categories 
allows us to directly compare their average SF activity, 
AGN probabilities, and perform X-ray stacking. Beyond 
this binary classification, we will also use quantitative es- 
timates of the Clumpiness, so as to distinguish the most 
clumpy galaxies from more moderate cases. 

The parent sample was selected in th e GOODS-South 
field, w ith deep HST/ACS imaging from lGiavalisco et al.l 
(|2004| ). This is the most suitable field combining deep- 
enough imaging for substructure identification (clumps), 
a large coverage of moderate-mass galaxies in optical 
spectroscopic surveys, with large-enough statistics. 

The main AGN identification tool used in this study 
is the MEx diagnostic. This narrow emission line diag- 
nostic is less affected by nuclear obscuration than X-ray 
selections, and was indeed shown to complement X-ray 
AGN selection with the identification of X-ray weak AGN 
(Jll) - X-ray stacking (Section 5) will confirm that most 
of the MEx-identified AGN in our sample are too faint 
to be individually detected in X-rays but are detected in 
stacked data. As opposed to the BPT technique (Bald- 
win, Phillips & Tcrlevich, 1981), the MEx diagnostic 
diagram remains useable beyond z = 0.4 with optical 
spectra, and its robustness was tested up to z w 1 in 
Jll. A key requirement for our sample selection is then 
the availability of a reliable [Oiii]/H/3 flux ratio mea- 
surement. The Blue diagra m ([Oii]A3727/H/?ew versus 
[Oiii]/H/3, lLamareiflell2010[) will also be used whenever 
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TABLE 1 

Properties of the galaxies in our Clumpy and Stable disk samples. 







Object 




Mass and redshift 


dumpiness 


Type 


ID 


a (J2000) 


<5(J2000) 


^■spec 




visual 


automated 




SI 


03:32:11.38 


-27:42:06.5 


0.733 


10.1 


1.4 


1.6 




S2 


03:32:17.63 


-27:48:11.8 


0.735 


10.7 


1.1 


2.6 




S3 


03:32:19.68 


-27:50:23.6 


0.559 


10.9 


2.1 






SI 


03:32:23.17 


-27:55:41.2 


0.733 


10.2 


2.3 


2.4 




S5 


03:32:23.40 


-27:43:16.6 


0.616 


11.1 


1.6 


2.1 




S6 


03:32:27.83 


-27:43:37.5 


0.548 


9.8 


1.2 


1.6 


Stable disks 


S7 


03:32:54.51 


-27:47:03.5 


0.533 


10.7 


1.2 






SS 


03:32:08.14 


-27:47:12.3 


0.578 


9.7 


1.6 






S9 


03:32:20.69 


-27:51:42.1 


0.679 


10.9 


2.4 


1.6 




S10 


03:32:30.43 


-27:51:40.3 


0.760 


10.4 


1.1 


2.7 




Sll 


03:32:29.17 


-27:48:33.1 


0.432 


10.1 


1.9 


1.4 




S12 


03:32:18.69 


-27:51:49.3 


0.457 


10.2 


1.4 


1.9 




S13 


03:32:19.78 


-27:54:09.1 


0.735 


10.4 


1.8 


3.1 




CI 


03:32:28.18 


-27:40:51.6 


0.426 


9.9 


1.3 






C2 


03:32:17.47 


-27:48:38.4 


0.737 


9 . 8 


4 i 


4.5 




C3 


03:32:19.61 


-27:48:31.0 


0.671 


10.3 


1.2 






C4 


03:32:36.68 


-27:39:54.6 


0.455 


10.2 


3.5 


3.1 




C5 


03:32:27.11 


-27:49:22.0 


0.559 


9.7 


1.6 






C6 


03:32:33.01 


-27:48:29.4 


0.664 


9.6 


1.8 




Clumpy disks 


C7 


03:32:34.04 


-27:50:09.7 


0.703 


10.2 


5.0 






C8 


03:32:51.52 


-27:47:58.1 


0.737 


11.0 


3.1 


4.7 




C9 


03:32:23.66 


-27:49:38.0 


0.578 


10.7 


1.1 


4.2 




CIO 


03:32:15.79 


-27:53:24.7 


0.676 


10.7 


1.7 


3.8 




Cll 


03:32:29.52 


-27:55:27.2 


0.663 


9.7 


1.2 






C12 


03:32:14.59 


-27:49:13.4 


0.562 


9.6 


1.6 


2.8 




C13 


03:32:15.35 


-27:45:07.0 


0.861 


10.5 


3.8 


2.7 




C14 


03:32:25.19 


-27:51:00.0 


0.841 


10.4 


4.0 


3.8 



TABLE 2 

AGN probabilities for our Clumpy and Stable disk samples, for various AGN categories, using the MEx and Blue 

DIAGNOSTICS SEPARATELY. 







Pmex 


^Bluc 


Type 


ID 


SF comp LINER Sy2 


SF comp LINER Sy2 





SI 





.49 





,27 


0. 


23 





.01 


0. 


17 





.82 


0, 


,00 





01 




S2 





.07 





19 





44 





00 





29 





,71 





00 





00 




S3 





23 





.74 





03 





00 





71 





.29 





00 





00 




SI 





.86 





.14 


0. 


00 





00 





.80 





.20 





00 


0. 


00 




S5 





.78 





22 


0. 


00 





00 





93 





.07 





00 


0. 


00 




S6 





92 





.05 





03 





00 





,44 





.26 





,27 





02 


Stable disks 


S7 





.15 





.71 





11 





00 





,68 





.31 





01 





00 




S8 





93 





02 





02 





03 





,58 





12 





02 





29 




S9 





03 





.40 





57 





01 





30 





.34 





34 





01 




S10 





12 





.40 





47 





01 





50 





19 





01 


0. 


00 




Sll 





43 





.31 





26 





00 





23 





31 





,44 





02 




S12 





.47 





.34 


0. 


19 





00 





,68 





29 





02 


0. 


00 




S13 





.89 





11 


0. 


00 





00 





91 





09 





00 


0. 


00 




CI 





04 





.01 


0. 


16 





80 





07 





.01 





00 





91 




C2 





80 





.03 





06 





11 





,57 





.09 





01 





32 




C3 





03 





.10 





53 





34 





43 





17 





00 


0. 


40 




C4 





11 





12 





47 





30 




















C5 





.85 





02 





03 





09 





43 





.03 





01 





53 




C6 





96 





.01 





01 





02 


















Clumpy disks 


C7 





02 





.05 


0. 


18 





15 




















C8 





00 





00 


0. 


21 





79 





00 





00 





,04 


0, 


95 




C9 





11 





.69 


0. 


20 





00 





,68 





30 





02 





00 




C10 





00 





.04 





70 





26 





47 





.14 





01 





38 




Cll 





43 





.04 





11 





42 





13 





00 





00 


0, 


87 




C12 





40 





.08 





04 





48 





01 





00 





00 





99 




C13 





09 





11 





16 





00 





32 





,67 





00 





00 




C14 





00 





.00 





16 





.84 





00 





,00 





04 





96 
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the ([Oii]/H/3)ew ratio can also be accurately measured, 
and the results from both AGN diagnostics will be com- 
pared. 

We thus select galaxies with optical spectra covering at 
least the H/3 and [Om] lines. We restrict this selection to 
the two largest optical spectroscopic surv eys in GOODS- 
South, namely GOODS /FO RS2 (IVanzella et all 
2001 and GOOD S/VIMOS dPopesso et all 120091: 
Balestra et all 120101) . For the latter, only spectra 
obtained with a medium-resolution VIMOS setup were 
considered because low-resolution setups would not 
yield accurate-enough flux measurements and in most 
cases cover too short wavelengths. The spectra used in 
this work hence all have a spectral resolution TZ > 720. 
We did not include other smaller spectroscopic datasct 
available in GOODS-South, which were examined^ but 
would add only 2 objects each to a sample of 27 galaxies 
while reducing the uniformity of the spectral coverage 
and offering lower signal-to-noise ratio spectra. 

2.2. Sample selection 

Among all objects covered in the GOODS/FORS2 and 
GOODS/VIMOS surveys, we considered only those com- 
plying with the following criteria: 

1. availability of a robust spectroscopic rcdshift z spe c 
with the hi ghest quality fl ags in the published 
databases (IVanzella et all 120081 : iPopesso et all 
[200l iBalestra et alll2010f ). 

2. redshift in the 0.4 < z spec < 1.0 range. The up- 
per limit is imposed by the availability of an [Om] 
detection, and also ensures that ACS imaging cov- 
ers optical rest-frame emission. The lower limit 
aims at obtaining two samples of clumpy and sta- 
ble galaxies that are about redshift-matchcd, con- 
sidering that clumpy galaxies are almost absent at 
z<0A. 

3. spectral range covering both the [Om] and H/3 
emission lines without being hampered by sky lines, 
a detection of both lines with a signal-to-noise ratio 
greater than four. 

4. stellar mass in the 9.6 < log(M»/M Q ) < 11.2 range: 
this arbitrary choice gives the largest possible sam- 
ple without introducing a strong mass bias between 
clumpy and smoother disks - we noted that in the 
studied redshift range, systems below our lower- 
mass limit are mostly clumpy, systems above our 
higher-mass limit are rarely clumpy. The technique 
employed for stellar mass estimates is described in 
Appendix [X] 

From this initial selection of 48 objects, we remove 
systems that cannot be obviously classified as Clumpy 
galaxies or Stable disks, namely: 

1. compact systems with major axis smaller than < 
0.8" in ACS i band, as no substructure could be 
distinguished (6 objects). 

3 This was done using the entire ESO 

CDFS spectroscopic compilation available at 

http : //archive . eso . org/wdb/wdb/vo/goods_CDFS_master 



2. interacting pairs and mergers with long tidal tails 
and/or double nuclei (4 objects). 

3. systems with major dust lanes hampering morpho- 
logical classification in optical bands (3 objects). 

4. spheroid-dominated systems, selected as having an 
axis ratio lower than 3 in ACS z and i bands with- 
out having clumps or spiral arms tracing a face-on 
disk (4 objects). 

5. systems harboring weak and/or elongated clumps, 
or clumps smaller than the ACS PSF, making the 
distinction between real clumps and short "floc- 
culent" spiral arms uncertain (4 objects). These 
weak clumps have a high contrast only in the bluest 
ACS bands (ultraviolet rest- frame), and merely vis- 
ible in the z band (optical rest-frame), which con- 
firms that they not similar to the massive clumps 
of Clumpy galaxies that are clearly visible in all 
ACS bands (see Fig©. 

Hence, 29% of the initial selection is rejected as these 
systems do not correspond to disk galaxies (compact ob- 
jects, spheroids, major mergers), and 14% of the sys- 
tems are rejected because they could not be classified as 
Clumpy or Stable disks (6%: major dust lanes, 8%:am- 
biguous or unresolved substructures). Representative ex- 
amples of rejected systems in each of the previous cat- 
egory are shown in Figure [5j This resulted in a parent 
sample of 27 systems with resolved morphology not cor- 
responding to major mergers, spheroids, and for which 
clumps/spiral arms are large/long enough to be distin- 
guished. We then classified these galaxies using the clas- 
sification criteria and dumpiness estimates detailed here- 
after, so as to build our final "Clumpy" and "Stable" 
disk samples (they are displayed on Figs. [3] and |4] after 
the classification is performed). 

Selecting galaxies with H/3 and [Om] detections can 
bias the sample toward AGN host galaxies. This would 
not hinder the comparison of Stable disks and Clumpy 
disks, but the absolute fractions of AGN in each sam- 
ple could be increased by this selection. However, the 
selection is mostly biased against passive galaxies. We 
aim at selecting star-forming disks, which should be on 
the Main Sequence of star formation, and the spectro- 
scopic surveys employed do typically detect these emis- 
sion lines in such galaxies and in the mass and redshift 
range that we consider, without requiring AGN. Indeed, 
even the objects with the lowest AGN probabilities in 
our sample, which re mostl likely star-forming galaxies, 
have [Om] and H/3 detected above our 4-tr threshold (and 
often above 5-cr), suggesting that similar star- forming 
galaxies would in general be selected. 

2.3. Morphological Classification and dumpiness 
Measurements 

In order to build two samples of "Clumpy" and "Sta- 
ble" disk galaxies, we performed dumpiness estimate us- 
ing an eyeball classification for the whole sample, and 
an automated measurement for about two thirds of our 
galaxies that have a robust disk/bulge luminosity model 
av ailable from the mo rphological study of GOODS-South 
bv lSalmi et all (l20H . 
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Visual classification — Examinin g the color images from 
the GOODS-HST/ACS images (IGiavalisco etalll200l . 
each of the authors independantly attributed a grade to 
each galaxy in the sample, between 1 and 5, using the 
following description of the grades: 

• 1.0: for robust smooth/stable disk, the stability 
being defined against bright clumps (spiral arms 
and bars, which are weaker instabilities typical for 
2 = disks, are not taken into account), 

• 2.0: for likely smooth/stable disk (some moderate 
clumps but spiral-dominated), 

• 3.0: for unsure systems, mix of stable and unstable 
regions, or potential low bulge/disk ratio, 

• 4.0: likely clumpy unstable disk (dominated by 
bright clumps but significant spirals), 

• 5.0: robust clumpy/unstable disk. 

In the following, we present results using the mean 
grade for each galaxy - median results were also exam- 
ined, and the mean grades were found to be slightly more 
conservative with respect to our final conclusions. We 
hereafter refer to the average grade for each galaxy as its 
visual dumpiness. 

Automated dumpiness measurement — A quantitative 
dumpiness measurement, based on bulge+disk axisym- 
metric luminosity model built with the Galfit software, 
is av ailable for a mass se lected sample in GOODS-South 
from iSalmi et al.l (|2012[ ) . Two thirds of our sample are 
covered by Salmi et al. - objects with a too low K- 
band luminosity or insufficient Galfit convergence are 
not covered. We then use these axisymmctric luminosity 
models, when available, to perform a dumpiness mea- 
surement which is slightly different from the one used 
initially by Salmi et al., in order to better distinguish 
clumps and spiral arms in the mass and redshift ranges 
that we are considering. 

For each galaxy, the ISalmi et al.l (|2012D model was 
built on the z-band ACS image within a segmentation 
map. The residual map is computed as the difference be- 
tween the image and model, and we divide the residual 
map by the model map to obtain the relative residuals. 
Negative residuals, which cannot correspond to bright 
clumps, but rather to inter-arm regions, badly subtracted 
bulges, or extended outer disks, are filtered out. We 
also filter out all relative residuals lower than 20%, be- 
cause they typically correspond to faint extended struc- 
tures: spiral arms or extended outer disks, not bright 
clumps. Such filtering was not applied in the Salmi et al. 
study. Here the 20% threshold value was found to give 
the best selection of "clumps" with respect to smoother 
structures, in particular spiral arms, as illustrated by the 
residual maps displayed in Figure Q] 

The automated dumpiness is defined by the sum of the 
residuals within the segmentation map divided by the 
number of pixels and multiplied by 1.8. We use the log- 
arithm of this value, so as to obtain values more directly 
comparable to visual estimates. The arbitrary multi- 
plication by 1.8 simply aims at obtaining similar scales 
from about 1.0 to 5.0 for both the visual and automated 
measurements. The correlations between dumpiness and 



AGN probability presented in the following Sections hold 
also without the use of a logarithmic scale. 

Selection of Clumpy disk and Stable disk samples — The au- 
tomated dumpiness measurements, when available, are 
tightly correlated with visual estimates (Fig. [5]). One 
single object shows a substantial difference (hereafter la- 
beled SI 3), caused by a strong spiral arm that was too 
strong to be filtered out from the relative residual map0: 
the visual estimate appears more relevant for this specific 
case. 

In order to obtain the largest possible samples of 
"Clumpy" and "Stable" galaxies, our main selection re- 
lies on the visual dumpiness estimates, the reliability 
of which was confirmed by the comparison with au- 
tomated measurements. We will nevertheless present 
results based on the automated measurement in Sec- 
tion 3.6. 

Here we classify the galaxies with visual dumpiness 
larger (respectively smaller) than 3.0 as "Clumpy" and 
"Stable" objects, which results in: 

• a sample of fourteen Clumpy disks displayed in 
Figure [3] and numbered CI to C14: their optical 
morphology is dominated by bright clumps, with 
or without spiral arms, similar to higher-redshift 
clumpy systems and suggesting a similar process of 
gravitational instability in a gas-rich disk. 

• a sample of thirteen Stable disks displayed in Fig- 
ure [4] and numbered from SI to S13: their opti- 
cal morphology is dominated by the usual struc- 
tures found in nearby disk galaxies, mostly bars 
and spiral arms. Star-forming clumps can be 
present, in particular along the arms, but are low- 
luminosity ones that do not dominate compared to 
the arm/inter- arm contrast, as in nearby late- type 
spirals. The presence of clumps only along spiral 
arms indeed indicates that such disks are overall 
stable, with Q ^ 1 instability reached only locally 
once the gas is compressed in the arms, as opposed 
to the very clumpy disks that are globally unstable 
with Q ;$ 1 instability not restricted to compressed 
spiral arm regions. 

We note that for objects with an automated dumpi- 
ness measurement, if we separate them using the median 
value of these measurements, they do not change cate- 
gory (see Fig. [2]). The only exception is galaxy S13: as 
already discussed, this specific case is dominated by a 
strong spiral arm and using its visual estimate is more 
relevant. 

Objects C6 and C7 w ere independently classified by 
lElmegreen et al.l ()2009al ). in agreement with our present 
classification. These authors also noted that the last ob- 
ject in our Figured] (rejected from the sample) does not 
show typical clumps, but rather low-mass structures that 
could be a combination of short flocculent spiral arms 
and dust absorption patches. Ob ject S3 was c l assifie d 
as a regular spiral disk galaxy bv iNeichel et al.l (|2008| ). 
consistent with our classification. 

4 unless a higher filtering threshold is employed, but this would 
eliminate clumps in several other objects. 
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Fig. 1. — Representative examples of ACS z-band images and associated relative residuals maps, using the Salmi ct al. (2012) luminosity 
profile models, after filtering of the relative residuals lower than 20%, showing that spiral arms, bars and bulges are mostly removed by this 
20% threshold, while bright clumps are kept, independent of being isolated clumps or clumps along spiral arms. The last object, labeled 
S13 in the following sections, is the only clear case of contamination by a spiral arm that was too strong to be filtered out. The galaxies 
shown here are labelled C13, Sll, C2, C9, C4, S10, S4 and S13, respectively, in the following sections and figures. 

The Stable disk sample contains 13 galaxies, at a median 
spectroscopic redshift of 0.62 (quartiles: 0.54 and 0.73) 
and with a median log(M*) = 10.42 (quartiles: 10.1 and 
10.8). Thus the two samples are relatively mass-matched 
and redshift-matched. Our samples consist of moderate- 
mass galaxies, although they also contain galaxies more 
massive than today's Milky Way. As previously detailed 
in Section 2, this is the expected redshift evolution of 
the violent instability process in galaxies fed by rapid 
gas inflow. 
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Fig. 2. — Correlation between the visual and automated dumpi- 
ness measurements. The same samples of Clumpy and Stable disks 
are obtained by using the median value of either measurement - 
one single object (hereafter S13) has a somewhat higher automated 
clumpiness value, known to be caused by a badly subtracted spiral 
arm, and is considered as a Stable disk. 

Some galaxies in the Stable disk sample have bright 
patches or moderate clumps along their spiral arms, such 
as S3 or S9 (as is also the case for many nearby spi- 
rals in near-UV or blue bands). When we use a binary 
Clumpy/Stable classification, such objects belong to the 
Stable disk sample, based on their clumpiness measure- 
ments. At some point we will also compare the presence 
of AGN with the individual clumpiness value for each 
galaxy, which will fully reflect the fact that some for 
galaxies of the "Stable" sample do have some moderate 
clumps. 

2.4. Main properties of the Clumpy disk and Stable disk 

samples 

Mass and redshift distribution — The Clumpy disk sample 
contains 14 galaxies, at a median spectroscopic redshift 
of 0.66 (quartiles: 0.56 and 0.74) and with a median 
stellar mass log(M*) = 10.2 (quartiles: 9.8 and 10.6). 



Star formation rates — Clumpy disk galaxies generally 
have relatively high star formation rate per unit stellar 
mass (sSFR, specific star formation rate), higher than 
stable disks and spheroid-dominated systems, but gener- 
ally not as high as starbursting systems for which sSFR 
>> 1 Gyr -1 - they lie within the "Main Sequence" 
of s tar formation defin ed by Elbaz et al. (2011, sec 
also iNordon et al.l 1201 If) that dominat es the budget of 
star formation (Rodigh iero et al1l2011[ ). Relatively high 
but not starbursting sSFRs were previously observ ed for 
z ~ 2 clumpy disks (jForster Schreiber et al.l 120111 ) , and 
our intermediate-redshift sample has the same property, 
as shown on Figure [7] 

Individual detections in the 16 and/ or 24 ^m Spitzer 
data from iTeplitz et al.l (|2011l) and iDickinson et al.l 
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2004 ). converted into sSFR using the I Chary fc Elbaa 
2001f ) relations, indicate higher sSFR in Clumpy types, 
but only about two thirds of our targets are individu- 
ally detected, with large uncertainties. To obtain a more 
representative estimate of the average sSFR in each sam- 
ple, we stacked the 16 and 24 (im data, using the same 
stacking methodology as for X-ray stacking, detailed in 
Section 5. None of our targets lies within the Spitzer PSF 
of another identified infrared source. To avoid contami- 
nation by nearby sources, we measured fluxes within the 
FWHM of the PSF before correcting for the full aperture. 

Using these near-infrared estimates, we obtain an av- 
erage sSFR in our Clumpy disk sample of « 0.38 Gyr" 1 , 
somewhat higher than in our Stable disk sample (sSFR 
w 0.26 Gyr -1 ). Such moderately-high sSFRs in Clumpy 
disks compared to Stable disks are also supported by 
SED fitting results (see Appendix) and by the [On] equiv- 
alent widths and H/3 fluxes for Clumpy and Stable disks 
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Fig. 3. — Sample of 2 ~ 0.7 Clumpy disks. Images are 6x6 arcsec, from HST/ACS B, V and i bands. Spiral arms may be present but 
as not as contrasted as the main clumps. 
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Fig. 4. — Same as Fig. [3] for our sample of Stable disks. Disk morphologies are dominated by spiral arms and bars, with weak clumps 
forming only in the spiral arms, indicating global stability with only local instabilities in the arms, as opposed to the "Clumpy disk" sample. 
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Fig. 5. — Examples of rejected systems in our classification, because of (1) major axis smaller than 0.8 arcsec, (2) evidence for merger 
(double nuclei and/or tidal tails), (3) major dust absorption, (4) low axis ratio without spiral arms, suggesting spheroid-dominated nature, 
(5) substructures of unclear or unresolved type with weak and/or elongated clumps (two examples shown): in the first example, dust lanes 
seem to cut across spiral arms, resulting in apparent structures that could be either elongated clumps or short arms, in the second case the 
substructures are not resolved by the ACS PSF, making their nature (clump vs. short arms) unclear. In such rejected cases, the putative 
clumps arc weak and visible only in the bluest ACS bands (i.e., ultraviolet rest-frame, see Fig.(6]l. 



Clumpy disks (clumps in all HST/ACS bands) Rejected (no/weak clumps in redder bands) 

1 * / • % • 
\% / • * • 

Fig. 6. — Comparison of systems classified as "Clumpy disks" and rejected cases in two HST/ACS bands: the V band (F606W, top) 
corresponding to near-UV emission for our targets, and the z band (F850LP, bottom) corresponding to optical emission around the V 
band rest-frame. Systems classified as "Clumpy disks" have major clumps in the optical rest-frame, even if slightly less contrasted than in 
the UV rest-frame. Rejected systems have weak elongated clumps that are strong in the UV rest-frame but merely visible in the optical 
emission. All images use the same logarithmic grey scale. 
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Fig. 7. — Average or median position of our two samples in 
the redshift-sSFR plane, for various estimates of the star forma- 
tion rate: SED fitting (median values for objects with SED fit- 
ting, see appendix), 24 fim detections converted into SFR using 
Chary & Elbaz ( 2001) (median values for about two third of each 
sample) , and mean values from 24 /jm stacking (average values 
for the whole samples). The dash-dotted lines show the limits of 
the "M ain Sequence" of star formation, as defined in iElbaz et alj 
(120111 ), the dotted line shows the average Main Sequence. Clumpy 
galaxies have higher sSFR then Stable disks, but remain within the 
Main Sequence, not in the "starburst" mode, as expected for gas- 
rich, gravitationally unstable disks. Note that only 24 fim stacking 
cover the entire samples, while 24 fim detections and SED fitting 
are available for parts of the samples, hence lying at somewhat 
different average/median redshifts. 

(after estimating and subtracting the AGN component, 
see Section 3.7) . We show in Figure[7]the position of our 
two samples in the redshift-sSFR plane for these various 
measurement techniques. This confirms that we have 
selected a sample of gas-rich disks with sSFRs within 
the upper part of the "Main Sequence" of star formation 
at z ~ 0.7 and no t at the level of starbursting mergers 
(jElbaz et al.l 120 111 ) , and a comparison sample of more 
stable disks with lower sSFR, in the lower part of the 
"Main Sequence". 

The average star formation rate surface density of our 
Clumpy disks corresponds to a gas to total baryonic mass 
fraction of 35% (using the optical major axis size for the 
disk diameter, and inverting the star formation - gas sur- 
face density relation for disks from Daddi et al. 2010b) 
. Similar estimates yield a typical gas fraction 23% for 
the Stable disk sample. Applying the stability calcu- 
lations presented in Bournaud & Elmegreen (2009), we 
determine that violent disk instability (i.e., a Toomre pa- 
rameter for gas and stars Q < 1 over most of the disk) 
requires a gas fraction 30%, for the average stellar 
mass and optical size of our sample galaxies, and assum- 
ing a bulge fraction of 20%. The gas fractions estimated 
above are consistent with this instability threshold. This 
confirms again that we have successfully selected a sam- 
ple of gas-richer, violently unstable disks (presumably 
representative of z^2 clumpy galaxies), and a sample of 
stable disks with lower gas fractions (more representative 
of low- redshift spiral galaxies) . 

3. EMISSION LINE DIAGNOSTICS 



In this section we use AGN diagnostics based on stel- 
lar mass and/or emission line ratios. The methods used 
for stellar mass and emission line ratio measurements are 
detailed in the Appendix. We first study the AGN prob- 
ability from the MEx diagnostic and Blue diagram, using 
our binary classification in Clumpy/Stable sample. We 
will next study the correlation between AGN probabil- 
ity and Clumpiness measurement for individual objects 
(§4-5). 

3.1. Optical spectra 

While we will use AGN diagnostics based on the indi- 
vidual spectra of each galaxy, we first show stacked opti- 
cal spectra for the samples of Clumpy Disks and Stable 
Disks, so as to highlight the general properties of these 
two samples. We kept only galaxies for which the spec- 
tral coverage encompasses [On] as well as [Oni] and H/3 , 
and which do not have the [Oni]A4959 line affected by a 
strong sky line. This led to rejecting Clumpy systems C2, 
C4 and C7, as well as Stable disks S4 and S7: 78% of the 
sources are included in the stacked spectra, and the ex- 
cluded systems do not have extreme values of line ratios 
and AGN probabilities for their morphological class (see 
Table 2). Each individual spectrum was re- normalized to 
the same average continuum level in the 3800-4800 A rest- 
frame spectral range. Stacked spectra were normalized 
to the same number of sources per sample. The stacked 
spectra, smoothed with a Gaussian kernel of FWHM 5 A , 
are shown on Figure |SJ 

The stacked spectra show that the H/3 emission, 
normalized to the same average continuum level, is 
marginally («20%) stronger in clumpy galaxies than in 
stable disks (consistent with a somewhat higher sSFR). 
The underlying absorption is weak in both cases, slightly 
higher in stable disks (as confirmed by H7 and H<5 lines). 
The [On] emission is stronger by a factor 1.8 in the 
stacked clumpy disk spectra, while the [Om] emission 
in clumpy types is a factor 5.2 stronger in clumpy disks, 
and the [Oni]/H/3 emission flux ratio higher by a factor 
3.8. The continuum shape is quite similar in both cases. 
In the following sections, we use the MEx diagram and 
other diagnostics to show that the strong [Oni] excitation 
by clumpy galaxies is a likely AGN signature, and cannot 
be explained simply by a low gas metallicity effect. 

3.2. MEx diagnostic 

Figure [9] shows our clumpy disk and stable disk sam- 
ples on the MEx diagram ([Om]/H/3 versus M», Jll). 
The [Om] excitation is clearly higher in clumpy disks, 
and these mostly lie above the empirical dividing line for 
AGN host galaxies, defined in Jll. Stable disk galaxies 
have lower [Om] excitations and arc mostly identified as 
non-AGN galaxies on the MEx diagram; only a few high- 
mass cases lie in the MEx-intcrmediate or MEx- AGN re- 
gions (respectively between and above the dividing lines 
on Figure H|) . 

In addition to the empirical dividing lines on Fig. [9j 
the MEx diagnostic was calibrated using the BPT clas- 
sification of >10 5 SDSS galaxies in order to quantify the 
probability to observe an AGN host system as a function 
of [Om]/H/3 and M* (see Jll). The SDSS calibration 
galaxies lie at low redshift (z < 0.1) and were first classi- 
fied using standard BPT diagrams into the following cat- 
egories: star- forming (SF), composite (comp), LINER, or 
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Fig. 8. — Stacked optical spectra for our samples of "Clumpy" and "Stable" disks, normalized to the same continuum level in the 3800- 
4800A spectral range and to the same number of galaxies per sample, and smoothed with a Gaussian kernel of FWHM 5A. Wavelengths 
at which significant sky lines contaminate a single spectrum, and thus potentially the stacked spectrum, are indicated: the main emission 
lines arc left unaffected. 
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Fig. 9. — Clumpy and Stable disk samples on the MEx diagram. 
The position of the stacked spectra from Figure \E\ is also shown 
(stars). The dividing lines separate the MEx-SF, MEx-AGN, and 
MEx-intermcdiate regions defined in Juneau ct al. (2011). The 
gray-scaled background colors indicate the AGN probability, cali- 
brated using the BPT diagnostic on z ~ 0.1 SDSS galaxies. The 
evolution of the mass-metallicity relation out to z R3 1 in the stud- 
ied mass range should not affect line ratios by more than 0.1 dex 
(see text Section 4.4). 

Seyfert 2 (Sy2). The likelihood of a certain spectral class 
is defined as the relative fraction of SDSS galaxies of that 
particular class within the lcr uncertainties on the MEx 
diagram. The four spectral classes listed above are mutu- 
ally exclusive so their sum is, by definition, equal to unity 
(P(SF)+P(comp)+P(LINER)+P(Sy2)=l). Individual 
AGN probabilities on the MEx diagnostic are indicated 
in Table [1] As discussed later in Section l3T4l the proba- 
bilities should not be strongly affected by the metallicity 
evolution between z < 0.1 (the rcdshift range calibration 
sample from Jll) and z ~ 0.7 (the typical redshift of our 
samples), especially in the mass range studied here. 

Stable disks have a median Pmex(SP) of 47% and, 
except for 2-3 AGN candidates, their AGN probabili- 
ties are low. In particular, the probability of hosting a 
Seyfert 2, Pmex(Sv2) is always below 3%, showing that 
systems that may not be pure star-forming galaxies show 
only weak signs of BH activity, with either composite- like 
properties or LINER-type excitation. 

The MEx probabilities for Clumpy disks are quite dif- 
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Fig. 10. — Bl ue diagram with the dividing lines defined in 
(Lamarcillc 2010) four our galaxy samples. The background prob- 
abilities are defined and coded as on Figure [9] The [On] /H/3 ratio 
is here an equivalent width ratio, while [On] /H/3 is a flux ratio. 

ferent. Except for three systems where Pmex(^4GA^) is 
low, the probability for the measured [Oni] excitation to 
result purely from star formation is always below 50%, 
and often close to 0. Half of the clumpy galaxies have a 
PMEx(<Sy2) comparable to or higher than the probability 
of composite properties or LINER-type excitation. 

3.3. Blue diagram 

We use the Blue diagram (|Lamar cillc 201(|) to perform 
a second classification of our samples. In some cases, the 
Blue diagram may be less efficient than the MEx diag- 
nostic in separating star-forming and AGN host galaxies 
because of an overlap region gathering galaxies identified 
as star-forming and as Sy2 on the BPT diagnostic. On 
the other hand, the Blue diagram may more clearly sep- 
arate LINER excitation (with high [On] excitation) from 
Sy2 activity (with the highest excitation in [Oni] ). 

Our two samples are shown on the Blue diagram on 
Figure [TUl and Table 1 gives the associated probabilities 
(calibrated on SDSS galaxies, as already explained for 
the MEx diagnostic). This confirms that a large fraction 
of clumpy disks likely host an AGN, unlike the control 
sample of stable disks. Furthermore, the figure shows 
that the driving difference is the [Oni] excitation, and 
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Fig. 11. — (R23;032) diagram (see text for definitions), showing 
the position of our two samples compared to the distribution of 
SDSS galaxies in four spectral classes: star-forming, composite, 
LINER and Seyfert 2, respectively, based on the BPT diagnostic 
applied to SDSS galaxies as in Jll. We did not apply any dust 
correction to our samples and to SDSS galaxies, as we estimate 
the correction to O32 to be lower than 0.1 dex. 

that clumpy disks are characterized by Sy2 or SF+Sy2 
activity, rather than LINER excitation. This feature is 
also apparent on the individual classification probabili- 
ties listed in Table 1, where the values of Psiue (LINER) 
are all very small. 

3.4. A metallicity effect? 

Our Clumpy galaxy sample shows, for most sources, 
the expected properties of AGN host systems (together 
with with star formation) on both the MEx and the Blue 
diagnostics. However, an alternative explanation for high 
[Onil/H/3 flux c ould be a low gas phase metallicity (e.g., 
iMcGaughl 119911 ). Theoretically, the disk instability as- 
sociated to clumpy morphologies is interpreted as being 
linked to the infall of fresh, low-metallicity gas (see in- 
troduction, DSC09 and Bournaud & Elmegreen 2009), 
so these systems could a priori be mctal-dcfficicnt. Nev- 
ertheless, simulations indicate that the gas is rapidly en- 
riched in the dense star-for ming clumps. Analysis of cos- 
mological simulations from iCeverino et al.1 ()2010P) gener- 
ally reveals solar metallicity in the clumps and marginally 
lower metallicities only in the inter-clump gas, typically 
one-third to half solar. Also, observations of clumpy 
galaxi es at z ~ 2 do not find them to be metal -deficient 
(e.g., iGenzel et al1[2008t iBournaud et al.ll2008h . 

As for our own sample of galaxies, the MEx diagnos- 
tic (Table 1) indicates that the probability to have the 
strong [Oin] excitation caused by a low metallicity ef- 
fect is, in general, quite lower than the AGN probabil- 
ity. These probabilities are calibrated on SDSS galaxies 
at z < 0.1 (Jll), however the redshift evolution of the 
mass-metallicity relation up to z ~ 0.7 is not sufficient 
to explain the high [Om]/H/3 ratio that we observe in 
Clumpy disks. For the mass range that we consider here, 
the typical variation of 12 + \og(0/H) from z < 0.1 to 



z r~j 0.7 is at most 0.1-0.2 dex, as shown by the Gemini 
Deep Deep Survey at 0.4 < z < 1.0 by iSavaglio et al.l 
(|2005l) . These authors found a slow evolution of the 
mass-metallicity relation with redshi ft at z < 1, ex- 
cept for lower- mass galaxies (see also iCresci et al.l 120121 
on the limited evolution of the mass-metallicity below 
redshift one). The corresponding variation of R23 = 
([Oiii]AA4959,5007+[Oii]AA3726,3729/H / 3 is not larger 
than 0.1 de x (using fo r insta nce the R23-metallicity re- 
lations from iMcGau gh (199^0) at fixed ionization param- 
eter). This possible redshift evolution is insufficient to 
strongly affect the position of our samples on the MEx 
diagram, nor the computed AGN probablities, and is ac- 
tually well within the assumed uncertainties. The com- 
parison to the Stable disk sample at similar redshift con- 
firms that the high [Om]/H/3 ratios do not result from 
redshift evolution in the mass-metallicity relation in the 
studied mass range, as such an effect would affect both 
samples. 

Furthermore, Clumpy galaxies are different from 
purely star-forming galaxies of any metallicity in the 
SDSS. This is shown by the (R 2 3;0 32 ) diagrams (Fig. El, 
where R 23 =([Oiii] AA4959, 5007+ [On] AA3726, 3729)/H,g, 
and 032=[Oiii]/[Oii]. The position of our Clumpy and 
Stable disk samples is compared to SDSS galaxies iden- 
tified as star-forming, composite, LINER and Seyfert 2 
on the BPT diagnostic (following the analysis performed 
in Jll). The high R23 values for Clumpy galaxies come 
along high O32 ratios. Note that high values of R23 do 
not trace low gas-phase metallicities in case of AGN exci- 
tation. Most Clumpy disks occupy the region of the pa- 
rameter space best represented by Seyfert 2's, and do not 
have the properties expected for star-forming galaxies of 
any metallicity, or for LINER excitation. Here again, 
the redshift evolution of the mass-metallicity relation in 
the studied mass range cannot explain this effect. Purely 
star-forming galaxies in our mass and redshift range lie 
almost exclusively in the R23<5 and 032<1 area of the 
(R 2 3;Q32) diagram, whatever their metallicity (as also 
shown bv ISavaglio et aLll2005h . The majority of Clumpy 
disks have higher R23 and O32 ratios, and are hence dif- 
ferent from purely star-forming galaxies of any metallic- 
ity. 

3.5. [NeIII]\3869 excitation 

A few other AGN diagnostics based on o ptical lines 
have been proposed in the literature (e.g. lYan et afl 
120111 ). A particularly useful one to complement the 
MEx diagnostic, given the spectral coverage available 
for our sample, is the recently proposed TBT diagnos- 
tic (jTrouille et al.ll201lD that compares the [Nem] A3869 
/[On] excitation to the (g — z) rest-frame color. This 
diagnostic typically requires to probe [Nem] /[On] ratio 
of the order of 0.1 or greater: the spectral data available 
for our sample are not sensitive enough for this in indi- 
vidual targets, but the [Nem] line is robustly detected in 
stacked spectral data. 

The TBT diagram for our Stable Disk and Clumpy 
Disk samples is shown on Figure El we find that Sta- 
ble Disks are dominated by purely star-forming systems, 
while a larger fraction of Clumpy Disks (but not necessar- 
ily all of them) should contain AGN, resulting in average 
pr operties on the AGN side of the AGN/SF limit defined 
bv lTrouille et al.l (|201lD . It is expected from the MEx di- 
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Fig. 12. — AGN probability versus Clumpiness value for Composite types (left), LINER types (middle), and Seyfert 2 types (right), using 
visual clumpiness values (circles) and automated clumpiness measurements (triangles). AGN probabilities in this figure are based on the 
MEx diagnostic. We indicate the corresponding linear correlation coefficients, and those based on the Blue diagram, in Table [3] 

Wc now consider the quantitative correlation between 
AGN probability and Clumpiness for all individual ob- 
jects, independent of the classification in two sub-samples 
of Clumpy and Stable disks used previously. 

We compare on Figure IT21 the AGN probabilities of our 
objects to their visual and automated Clumpiness esti- 
mates. This is done for the three considered AGN types 
separately (Composite, LINER and Seyfert 2). This ex- 
periment was also performed for combinations of AGN 
types (such as LINER+Seyfert 2, etc), considering either 
the MEx diagnostic or the Blue diagram, and either the 
visual or the automated clumpiness estimates. In each 
case we determined the linear correlation coefficients a, 
b and r for: 



TBT-SF 



■1.5 



-1.0 -0.5 
log [Nelll] / [Oil] 



Fig. 13. — TBT diagram for Clumpy and Stable disk samples 
(mean stacking) and for similar sub-samples weighted by the MEx- 
based probability to host AGN (for Clumpy disks) or to be purely 
star- forming (for Stable disks). 

agnostic that not all Clumpy disks have optical line exci- 
tation dominated by AGN, and hence that their average 
properties are intermediate between purely star-forming 
galaxies and AGN-dominated ones, as found here on the 
TBT diagram. 

To further confirm the presence of AGN in many 
Clumpy galaxies, we have studied the stacked data for 
the Clumpy Disk sample with statistical weights pro- 
portional to Pme*(AGN), the MEx-based probability of 
containing AGN, for each object. The idea is equiva- 
lent to selecting only the objects most likely to contain 
AGN among the Clumpy sample, improved by the use of 
statistical weights. The results are shown on the TBT 
diagram (Fig. fl"3|) . and the weighted data arc clearly in 
the TBT-AGN region. Hence, while the whole sample 
of Clumpy Disks has intermediate properties suggesting 
a mix of star-forming and active galaxies, the Clumpy 
disks identified as likely AGN hosts by the MEx diagram 
are confirmed to be dominated by AGN with the TBT di- 
agram. Conversely, the stacked data for the Stable Disk 
sample weighted by Pme*(SF) = 1 — Pmex{AGN), i.e. a 
weighted selection of Stable Disks supposed to be likely 
purely star-forming according to the MEx, lie clearly in 
the TBT-SF region on Figure [H 

3.6. Quantitative correlation between AGN probability 
and Clumpiness 



P(AGN) = axC + b ±r 



(1) 



where P(AGN) is the AGN probability for the consid- 
ered AGN type, C is the clumpiness value, and r is the 
Pearson correlation coefficient. Results are given in Ta- 
ble El 

The global correlation between AGN probability and 
clumpiness is dominated by Seyfert 2-typc activity: the 
Seyfert 2 probability vs. clumpiness correlation is the 
tightest one and drives the global AGN probability vs. 
clumpiness correlation, much more than LINER and 
composite types. This result is independent of the diag- 
nostic used (MEx or Blue diagram) and preferred clumpi- 
ness measurement (visual or automated). The most 
clumpy objects have spectral properties that are weakly 
more typical of Composite galaxies compared to the less 
clumpy ones, and not more typical (or even less typical) 
of LINER galaxies, but they clearly have spectral prop- 
erties that are more typical of Seyfert 2-like activity. 

These results confirm the previous conclusions, inde- 
pendent of the chosen clumpiness threshold used to sep- 
arate the Clumpy and Stable sub-samples, and of the 
reliability of visual clumpiness estimates. Furthermore, 
the fact that the strongest correlation is obtained when 
Seyfert 2-type activity alone is considered, and the ab- 
sence of a significant correlation between LINER activ- 
ity and clumpiness, confirms that the observed activity 
in clumpy galaxies corresponds to BH accretion, rather 
than shock-induced excitation or low-metallicity star for- 
mation that might have resulted in LINER-like spectral 
properties - which was unlikely, as supported by the pre- 
vious discussion on the evolution of the mass-mctallicity 
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TABLE 3 

Linear correlation coefficients between visual/ automated clumpiness values and AGN probability, for various 

COMBINATIONS OF AGN TYPES. 



AGN type(s) 




MEx 






Blue 




a 


b 


r 


a 


b r 



Using visual clumpiness: 














Composite + LINER - 


- Seyfert 2 


5.64 


45.3 


0.05 


8.31 


33.6 


0.16 


Composite - 


- Seyfert 2 


3.46 


29.6 


0.04 


10.4 


22.6 


0.37 


LINER - 


- Seyfert 2 


13.0 


1.90 


0.26 


16.8 


-14.4 


0.39 




LINER 


2.16 


15.8 


0.02 


-2.07 


11.0 


0.06 




Seyfert 2 


10.81 


-13.9 


0.31 


18.9 


-25.4 


0.47 


Using automated dumpiness: 














Composite + LINER - 


- Seyfert 2 


8.21 


54.9 


0.55 


2.19 


46.8 


0.07 


Composite - 


- Seyfert 2 


9.35 


11.8 


0.17 


8.30 


26.3 


0.28 


LINER - 


- Seyfert 2 


14.49 


-0.9 


0.18 


14.20 


-11.9 


0.39 




LINER 


-1.13 


28.2 


0.01 


-6.11 


24.5 


0.45 




Seyfert 2 


15.63 


-29.07 


0.57 


15.49 


-29.4 


0.57 



relation. 

It appears on Figure [12] that low-clumpiness galaxies 
have Composite-like properties rather than Seyfert 2- 
like properties, suggesting that if they contain AGN, 
these are lower-luminosity ones that do not dominate 
the [Oiirj/H/3 ratio compared to star formation, as op- 
posed to Clumpy disks where the [Oin]/H/3 ratio can 
be AGN-dominated, leading to high Seyfert 2-like inte- 
grated probabilities. 

3.7. Summary: optical line properties of Clumpy AGN 
candidates 

Clumpy galaxies show a much higher [Oin] excitation 
than Stable disk galaxies in the same mass and redshift 
range. Comparisons with the properties of [On] and 
[Nein] indicate that this trend is unlikely to be caused 
by mctallicity effects in purely star-forming galaxies, and 
that the AGN fraction should be higher in Clumpy galax- 
ies instead. Using the MEx diagnostic, not all Clumpy 
galaxies need to contain AGN, but many of them have 
a high AGN probability and the AGN probability corre- 
lated with both visual and quantitative estimates of the 
clumpiness. Similarly, some Stable disks may contain 
AGN, but the estimated fraction is significantly lower. 

We list the average [Oin] , [On] and H/3 luminosities for 
both samples in Table |4] Overall, Clumpy disks have an 
H/3 luminosity higher by about 20% than Stable disks, 
but their [Oin] luminosity is higher by a factor five. This 
is not directly representative for the properties of the 
hosted AGN because some objects in the Clumpy sam- 
ple are unlikely to host AGN, and vice versa. To better 
probe the properties of active clumpy disks, we have de- 
fined a "Clumpy AGN" sample, for which average prop- 
erties are measured over our Clumpy disk sample with 
statistical weights equal to Pmcx(AGN), i.e. we give the 
highest weights to the systems that are the most likely to 
contain AGN. Conversely, we defined a "Stable SF disk" 
sample, which corresponds to the Stable disk sample with 
statistical weights proportional to the MEx probability of 
being purely star-forming. This allows a more relevant 
estimate of the properties of AGN hosted by Clumpy 
disks (see Table [4j. In particular, the typical [Oin] lumi- 
nosity of Clumpy AGNs is 2.1 ± 0.2 x 10 41 ergs" 1 . This 
is ten times larger than in Stable SF disks, a result con- 
sistent with the more general properti es of AGN hosts 
and star-forming galaxies in the SDSS (jkauffmann et all 



120031) , and suggesting that ~ 90% of the [Oin] luminosity 
in Clumpy AGN comes from the AGN itself. 

Hence, the average AGN luminosity in the clumpy 
AGN candidates is L[om],AGN ~ 1-9 ± 0.2 x 10 41 erg s _1 

- corrected for the minor contamination by star forma- 
tion, but not corrected for dust extinction. This sug- 
gests an average bolometric AGN luminosity iboi.AGN ~ 
2 x 10 43 erg s" 1 (e.g. lNetzerll2009( ) for the Clumpy disks 
with high AGN probabilities. 

The average H/3 luminosity is slightly higher in the 
Clumpy AGN sample compared to the Stable SF one, 
by about 25%. Given that Clumpy disks have slightly 
lower SFRs (Section 2.4), this indicates that about one 
third of the H/3 luminosity in Clumpy AGN systems is 
powered by the AGN. The powering of the H/3 luminosity 
by star formation is thus about similar in Clumpy and 
Stable disks, with an extra contribution to the H/3 flux 
in Clumpy AGNs. 

4. X-RAY STACKING 

4.1. X-ray properties of individual sources and stacking 

AGN candidates in unstable disk galaxies are expected 
to have modest X-ray luminosities and significant obscu- 
ration by the high ISM column densities. They could 
thus be generally to faint to be individually detected 
even in deep Chadra surveys, especially in the hard band 
(Bll). Indeed, only one object in our sample, C10, is a 
known hard X-r ay selected AGN ([Alexander et al.ll2003l : 
iXue et aD 120 111 ). Several other sources are detected in 
the soft X-ray band of Chandra data (0.5-2keV), but in 
general with a signal-to-noise ratio that is too low to 
firmly distinguish between emission from an AGN or SF, 
both the soft X-ray luminosity and the star formation 
rate being relatively uncertain in these individual objects 

- object CI has a robust soft X-ray excess compared to 
its SFR, according to the Ranalli et al. (2003) relation, 
though. 

To confirm the high AGN frequency in clumpy galax- 
ies, we perfo r med mea n X-ray stacking as in, e.g., 
iWorslev et all (pOOol ) and lDaddi et~aTI ([2001 . All of our 
targets lie within 8' from the Chandra aim point in the 
CDFS, 80% of the sources being within 6.5', and the me- 
dian distance from the aim point being 4.3', where the 
sensitivity is relatively homogeneous (relative variations 
below 25%). No target lies within the Chandra PSF of a 
close (in projection) X-ray source. For each source, the 
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TABLE 4 

Optical line, X-ray and infrared luminosities for various sub-samples. We indicate mean values or values from mean 
stacking. Optical lines luminosities are uncorrected for dust extinction. Luminosities are in ergs -1 , except L ir in Lq. 



Mean luminosities 



MOIII] 
L [OIl] 

Lx, 0.5-8kcV rest. 
Lx, 2-10keV rest. 
^IR 



Stable SF disks* 

2.1 ±0.5 x 10 4 ~ 
3.3 ±0.4 x 10 40 

3.2 ±0.8 x 10 40 
1.0 ±0.3 x 10 41 
0.6 ±0.4 x 10 41 
5.0 ±0.5 x 10 10 



All Stable disks^ 
2.9 ± 0.4 x 10 4 ~ 
3.5 ±0.4 x 10 40 
3.8 ± 1.1 x 10 40 
1.2 ±0.3 x 10 41 
0.7 ±0.4 x 10 41 
5.2 ±0.4 x 10 10 



All Clumpy disks^ Clumpy AGN disksj 



1.5 ±0.1 x 10 41 
4.3 ±0.6 x 10 40 
6.9 ± 1.2 x 10 40 
4.2 ±0.4 x 10 41 

2.6 ±0.8 x 10 41 
4.5 ±0.6 x 10 10 



2.1 ±0.2 x 10 41 
4.1 ±0.5 x 10 40 
7.3 ±1.1 x 10 40 
5.1 ±0.6 x 10 41 
3.0 ±0.8 x 10 41 
4.8 ±0.7 x 10 10 



a Stable disk sample with (1 — PyiEx(AGN)) weights applied to each object to compute average values and to stack X-ray and infrared 
data, hence statistically selecting objects that have the highest probabilities to be purely star-forming. 
b Entire Stable disk sample with equal weights for each object - including the few with high AGN probabilities. 
c Entire Clumpy disk sample with equal weights for each object - including the few with low AGN probabilities. 

d Clumpy disk sample with Pmex{AGN) weights applied to each object to compute average values and to stack X-ray and infrared data, 
hence statistically selecting objects that have the highest probabilities to be clumpy AGN hosts. 



X-ray data were multiplied by a factor proportional to 
the square of the luminosity distance - stacked results 
without this correction factor were also examined and 
yield similar results. The stacked data, presented in Fig- 
ure I14[ were also normalized to the number of sources 
per sample. 

4.2. Stacked X-ray detections 

In the soft X-ray band (0.5-2 keV), X-ray stacking 
yields a 4.5-er detection^] for the Stable disks sample, and 
a 3.7 ± 0.2 times larger signal for the sample of Clumpy 
galaxies. The stacked signal is not dominated by the few 
most luminous systems in soft X-rays: the results are 
similar if we select randomly only half of the galaxies in 
each sample. Hence the factor of sa 4 soft X-ray excess 
is a general property of our whole Clumpy galaxy sam- 
ple, not just of a few outliers within this sample. In the 
hard band (2-8 keV), there is no detection for Stable 
disks, and only a marginal detection for Clumpy disks. 
The stacked data show a 3.5-cr detection but a single ob- 
ject (C10, which is the individually detected source) con- 
tributes 35% of stacked the signal and rest of the sample 
gives only a 2.4-ct signal. 

There is thus a higher X-ray luminosity in Clumpy 
galaxies than in Stable disks. The hard/soft count ratio 
for Clumpy disks is H/S ~ 6.5 if the stacked signal in 
the hard band is robust - conservatively, this is rather 
an upper limit. In the following we assume a spectral 
slope r = 1.2 ± 0.5 for luminosity conversions, consis- 
tent with this (uncertain) H/S count ratio and typical 
for AGN that are relatively obscured but not Compton 
Thick (e.g., Alexander et al. 2003, Jll). The comparison 
of [Oni] and X-ray luminosities in our sample is indeed 
consistent with such substantial obscuration (see Sec- 
tion 5.4). The stacking results indicate an average X-ray 
luminosity L x(2 _ wkeVrest ^ = 2.6 ±0.8 x 10 41 ergs" 1 for 
Clumpy galaxies, versus 0.7 ±0.4 x 10 41 ergs _1 for Stable 
disks, both values being uncorrected for absoprtion (see 
Tabeirj}. 



4.3. Star formation rate and X-ray excess 

The X-ray luminosity of Clumpy galaxies is almost four 
times higher than that of Stable disks, but their star for- 
mation rates are not higher. They have somewhat higher 

5 (7 being measured from the background pixel-per-pixel fluctu- 
ations 



sSFRs, but lower stellar masses too, and their absolute 
SFR was found to be slightly lower than that of Stable 
disks (from SED fitting and 24 /xm fluxes, Section 3.4, 
also consistent with H/3 luminosities once corrected for 
AGN contribution, Section 3.7). 

SFR estimates from individual detections at 16 and/or 
24 /tm are on average 7.5M Q yr _1 for Stable disks, corre- 
sponding to Lx/SFR ~ 9 x 10 39 ergs _1 /M Q yr _1 , sug- 
gesting X-rays in these systems trace mostly star for- 
mation (jRanalli et all 120031 : [Grimm et al.l 120031) . In de- 
tail, this Lx I SFR ratio is somewhat above the main 
trend defined by Ranalli et al. 2003, suggesting that 
there might be weak AGN in some of the Stable disks, 
consistent with the MEx results. The average SFR 
estimate in Clumpy disks, at 6.1M Q yr _1 , results in 
Lx/SFR ~ 4 x lO 4o ergs-yM yr" 1 , indicating an X- 
ray excess in these systems, by a factor > 5 compared to 
the X-ray luminosity expected from their star formation 
activity using the Ranalli et al. relation. 

However, not all targets are individually detected at 
16 or 24 /tm. Therefore we stacked the 16 and 24 /mi 
Spitzer data for our two samples of galaxies. Here we ap- 
plied the same weighting with luminosity distance as for 
the X-ray stacking, but also checking that un-weightcd 
stacking gives similar results, and normalizing stacked 
data to the same number of object per sample. The re- 
sults are shown on Figure 1141 note that some infrared 
sources relatively close to our targets do appear on the 
stacked images, especially at 24 /im, but none of our tar- 
gets lies within the Spitzer PSF of another identified in- 
frared source, and flux ratios between our two samples 
were measured within the FWHM of the PSF to avoid 
contamination. 

The results clearly show that there is no 24 or 16 /im 
excess in Clumpy galaxies. Instead, their stacked signal 
is 32±4% lower at 16 /im and 23±3% lower at 24 /im than 
that of Stable disk galaxies. 24 /tm are not proportional 
to star formation rates in a rcdshift-indcpcndcnt way. 
Here the main bias between the two samples would result 
from the highest-rcdshift objects being Clumpy ones at 
z ~ 0.85, where PAH lines are redshiftcd in the Spitzer 
24 /im band and increase the observed flux at fixed SFR. 
This may explain why the Clumpy-to-Stable flux ratio is 
somewhat higher (smaller than unity but closer to unity) 
at 24 /tm than at 16 /tm - a small effect that could only 
strengthen our conclusion. 

Thus, star formation cannot explain the X-ray excess 
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Fig. 14. — Stacked data for the Clumpy and Stable disk samples, for soft X-rays, hard X-rays, 24/i m and 16/x m data, respectively. The 
stacked data were re-normalized to the same number of objects per sample, to the same luminosity distance for each individual object (see 
text for details), and background-subtracted. The same linear color bar is used for each pair of image, normalized to 1.0 for the brightest 
central pixel value for soft X-rays, 16 and 24 /im, and to 0.6 for hard X-rays. Note that the two bright areas at a ~ 0.4 — 0.5 sec in the 
soft X-ray stack for Stable disks are detections of nearby sources and not indicative of the typical noise level. The same applies in a more 
evident way to the 16 and 24 fim stacks with some nearby sources detected. 

ties is subject to large uncertainties given the potential 
effects of obscuration. 



in Clumpy disk galaxies. One could still envision that 
Clumpy galaxies have an X-ray excess for some reason 
else than hosting AGN, for ins tance shocks in the IS1V0. 
However, IRanalli et al.l (|2012D found that the local X- 
ray relations for star-forming galaxies hold at increasing 
rcdshift, where clumpy types should be increasingly fre- 
quent. 

To further test the consistency of the AGN interpre- 
tation of the X-ray excess with the previous optical line 
results, we performed the X-ray stacking measurements 
for the Clumpy galaxy sample with statistical weights 
proportional to PMex(AGN) (previously defined as the 
"Clumpy AGN" weighted sample) and the correspond- 
ing "Stable SF" disk sample with weights proportional 
to (1 - Pmcx(AGN)). As indicated in Table 0] the aver- 
age X-ray luminosity increases in the weighted Clumpy 
AGN sample, compared to the un-wcightcd Clumpy disk 
sample, while SFR estimates (Lin, Lnp, SED fitting) 
remain relatively unchanged. Hence, the X-ray excess 
in Clumpy disks is strongest in objects that have the 
highest probability to host AGN based on their nebular 
line properties, as expected if the X-ray excess is indeed 
caused by frequent AGN in Clumpy disks. This further 
demonstrates the agreement of the two methods (X-ray 
stacking and optical lines) toward the AGN interpreta- 
tion. 

4.4. Inferred AGN and BH properties 

The identification of AGNs in Clumpy galaxies with 
various techniques indicates that the AGN dominates 
the [Oin] line emission and X-ray emission, on average, 
in this galaxy sample. Here, we discuss the typical BH 
mass accretion rate from the observed AGN luminosity 
although the interpretation of [Oin] and X-ray luminosi- 

6 There is probably strong supersonic turbulence and shocks in 
clumpy disks, however, spiral arms in Stable disks typically corre- 
spond to strong shocks in interstellar gas, too. 



Luminosity and obscuration — The average X-ray lumi- 
nosity of the Clumpy AGN sample (weighted using the 
MEx probabilities to be as representative of actual AGN 
hosts as possible) is Lx, 2-iokcVrost. ~ 3 x 10 41 ergs~ 1 , 
after estimating and subtracting t he contribution of sta r 
formation with the relations from IRanalli et al.l (|2003|) . 
Direct conversion into AGN bolometric luminosity from 
the X-ray data alone is uncertain as the obscuration can- 
not be measured using a hard band to soft band ratio for 
our sample. 

The comparison with [Oin] luminosities in Table 0] 
shows that ix/^[oni] < 2, which means that these 
AGN could be substantially obscuredj] but not neces- 
sarily Compton thick. This indicates that the intrinsic 
X-ray luminosity of these AGN would be higher than 
observed, and wc infer from the observed [Oin] lumi- 
nosities an intrinsic X- ray luminosity in the range of 
Lx.int. ~ lQ^-^crgs- ^Netzer et al.ll2006l ). This is con- 
sistent with the models in Bll, predicting high gas col- 
umn densities (Njj ~ io 23 ~ 24 cm -2 from the general 
large-scale ISM, not counting extra obscuration from a 
potential AGN torus) on most lines of sight but reach- 
ing Compton thickness on a fraction of these. This is 
also consistent with the observation by Mullaney et al. 
(2011) of numerous X-ray AGN on the Main Sequence 
of star formation that could be obscured without being 
Compton thick. 

Bolometric luminosity and BH accretion rate — Estimating 
the AGN bolometric luminosity from the observed [Oin] 
luminosity is quite uncertain. For instance, the conver- 
sion factor depends strongly on the X-ray luminosity cor- 
rected for absorption (e.g., Netzer et al. 2006), which is 

7 because unobscured AGNs would have £x/^[Olll] ^ 10 
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uncertain here, and the £boi/£[om] values reported in 
the literature span a range from a few tens to several 
thousands. If we adopt a relatively conservative ratio 
for MEx-selected AGN of «300, the average AGN bolo- 
metric luminosity is of the order of a few 10 ergs -1 . 
Some studies have suggested higher conv ersion factors 
for similar [Quil luminosities, up to »3000 (lNetzerj|2009t 
iHeckman et"aLlT20p . 120051: iShen et all |201j, in which 
case the AGN bolometric luminosity may reach a few 
lO^ergs- 1 . 

If the mass to luminosity conversion efficiency is 10%, 
the AGN luminosity is Lboi ~ 10%?tic 2 , then the av- 
erage BH accretion rate is m ~ 10~ 2 — 10 -1 M Q yr" 1 , 
depending on the adopted conversion factors. For a star 
formation rate of a few M yr -1 , this corresponds to 
an rh/SFR ratio around 10 -3 or somewhat larger, as 
found more g enerally for BH growth in gas-rich galax - 
ies at z > 1 (jDaddi et al.l 120071 : IMullanev et"aLll2012b| ) 
and implying that these BHs co uld be growing toward 
the observed scalin g rela tions (jMagorrian et al.l 119981 : 
iMerritt & Ferraresell200ll ). If the violent disk instabil- 
ity process can be maintained over timescales of a few 
Gyr (or even less for the highest estimates of the bolo- 
metric luminosity) , the resulting BH mass could be a few 
10 7 M©, possibly representing most of the final BH mass 
in galaxies with stellar masses of a few 10 10 Mq. The 
average AGN luminosity remains sub-Eddington (even 
for the highest [Om]-to-bolomctric luminosity conversion 
factors), but short accretion bursts with higher Edding- 
ton ratios could be possible, for instance if giant clumps 
migrate and coalesce centrally. 

5. DISCUSSION: COMPARISON TO HIGHER- REDSHIFT 
DATA AND MODELS 

Our sample of clumpy galaxies at z-~0.7 has proper- 
ties that are globally representative of gas-rich unstable 
disks, which are increasingly common at higher rcdshift. 
As detailed in Section 2, their relatively high but Main 
Sequence star formation rate are indicative of disks with 
high gas fractions rather than mergers. This suggests 
that the results obtained for our intermediate-rcdshift 
sample are also representative of the numerous unstable 
clumpy galaxies at z > 1 , most of which are g as-rich sys- 
tems o n the Main Sequence. At redshift zss2. lGuo et al.l 
(|201lD observed in the Hubble Ultra Deep Field that 
some clumpy disk galaxies contain X-ray detected AGN. 
Having only a limited fraction of such AGN individually 
detectable in X-rays is fully consistent with our findings. 
Our narrow line and X-ray stacking results actually in- 
dicate that there should be more AGN of this type than 
those unveiled by i ndividual X-r a y det ections. In mas- 
sive disks at z « 3, ICresci et al.l (|2010f ) noted that the 
[Oiii]/H/3 ratio sometimes peaks in inner regions of the 
disk, either at the center of close to it, not in the outer 
disk as expected for the usual inward metallicity gradi- 
ents. In the light of our results, this might be another 
signature of AGN activity in clumpy rotating disks. 

The triggering of nuclear activity and BH growth by 
violent disk instabilities should then be a general pro- 
cess in stream-fed, gas-rich, unstable disk galaxies at any 
(higher) rcdshift, and in particular at the peak epoch 
of cosmic star formation activity around redshift two. 
Our estimates of AGN luminosities and BH accretion 
rates are consistent with the theoretical expectations 



from Bll, and over timescales of one to a few Gyrs this 
process could fuel a large part of the final BH mass in 
modcratc-mass galaxies, and potentially even in mas- 
sive objects. Hence, the violent clump instability in 
high-rcdshift galaxies is much more efficient in feeding 
a BH than the weaker bar instabilities frequent in low- 
redshift systems. However, bars are the most consid- 
ered instability mode in theoretical models to date (e.g., 
iFanidakis et al.l [201 If ) ■ Another difference is that high- 
rcdshift clump instabilities can grow a B H together with 
either a classical bulge or a pseudo-bulge (lNpguchi|[l999l : 
lElmegreen et al.l2008bHInoue fc Saitohll2012[ ). while bars 
could only grow pseudo-bulges. 

As disk galaxies fed by rapid cold gas accretion could 
remain unstable for a few Gyr, the duty cycle of BH 
growth could be high (>10%) in this phase. This does 
not require BHs to be permanently active in these galax- 
ies: some Clumpy disks in our sample have low AGN 
probabilities and/or no X-ray excess with respect to their 
SFR. As an illustration, if we assume that ~50% of Main 
Sequence star-forming galaxies at z=l-2 are clumpy disks 
(consistent with deep imaging and spectroscopic surveys, 
e.g. Elmcgrccn et al. 2007 and Shapiro ct al. 2008), 
and 50% of them have AGN with L x > 10 42 crgs -1 , 
this could contribute to an AGN fraction of roughly one 
fourth in Main Sequence star-forming galaxies (consider- 
ing AGNs with an intrinsic Lx>10 ergs -1 ). Our data 
and the models in Bll suggest that clumpy disks can host 
more luminous AGNs, maybe even at Eddington-limited 
levels, but with lower AGN fraction and duty cycles for 
such higher-luminosity populations. 

These findings could explain how AGN are fre- 
quently fed in galaxies that do not display the signa- 
tures expect ed for major mergers (IGrogin et al.l [20051: 
Gabor et all 120091: IKocevski et ail 1201 It IMullanev et atl 



2012al) . In this long duty cycle mode, with ubiquitous 



instabilities in high-redshift disks fed by rapid gas accre- 
tion, a large part of the BH mass in today's galaxies could 
have been assembled in AGN phases that have moder- 
ate luminosities and significant obscuration, most often 
with bolometric luminosities of ~ if) 43-44 ergs -1 . Short 
and bright phases could occur when dense clumps reach 
the nucleus and the AGN luminosity reaches QSO levels 
and the BH accretion rate increases, potentially reach- 
ing Eddington-limited levels (Bll). In this case, each 
given system would spend most of the time at moderate- 
luminosity AGN levels, but it would be possible that 
a larger fraction of the BH mass is assembled during 
shorter high- luminosity phases (or not). This is impossi- 
ble to study in our small sample, but is suggested by some 
observations (although lacking infor mation on the nature 
of host galaxies - e.g., ISoltanlll982D and some cosmolog- 
ical simulations (although n ot resolving the detailed ac- 
cretio n processes yet - e.g., iMartizzi. Tevssier k, Moorel 

HOED). 

The contribution of violent disk instabilities to the lu- 
minosity function of AGNs and mass function of BHs re- 
mains to be established. Here we have provided evidence 
that high-redshift disk instabilities trigger BH growth: 
this can be through BH formation in the giant clumps 
(Elmegreen et al. 2008), through giant clumps migra- 
tion bringing gas directly into a central BH, or through 
the more general gas inflow that giant clumps trigger in 
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high-rcdshift disks even when they don't reach the cen- 
ter themselves (Bournaud et al. 2011). This is similar to 
the known and observed triggering of BH growth by in- 
teractions and mergers (e.g., Ellison et al. 2011). But al- 
though mergers do trigger BH growth, selections of AGN 
host galaxies do not find a substantial excess of mergers 
with respect to control sample, implying that merger- 
induced process does not dominate the triggering of AGN 
in the Universe. The process of violent instability result- 
ing from rapid cold gas accretion seems to dominate the 
stellar mass assembly of galaxies compared to big merg- 
ers (Dekel et al. 2009, Brooks et al. 2009, L'Huillier et 
al. 2011), so this process could be more likely to sig- 
nificantly contribute to AGN triggering. Koccvski et al. 
(in preparation) find that clump instabilities are not a 
dominant trigger of AGN with observed X-ray luminosi- 
ties above 10 ergs" 1 . This remains consistent with our 
present findings, as our sample of fourteen Clumpy disks 
would include only one such AGN, with numerous AGN 
revealed by other means - narrow line diagnostics, X- 
ray stacking, and in a few cases a possible X-ray excess 
compared to the SFR but below 10 42 ergs" 1 . 

6. CONCLUSION 

The growth of Supermassive Black Holes (BHs) trig- 
gered by disk instabilities in high-redshift galaxies could 
be hard to observe directly at redshifts z > 1, especially in 
X-rays. According to numerical simulations, this mode 
would, most of the time, trigger moderate luminosity and 
obscured AGNs. In this paper, we have built a sample of 
unstable (Clumpy) disks and a control sample of Stable 
disks at intermediate redshift (z~0.7). This way, more 
reliable Narrow Line diagnostics can be used compared to 
higher redshift, and moderate X-ray luminosities are in 
easier reach of (stacked) observations. Our classification 
of Clumpy and Stable disks used visual estimates of the 
dumpiness and was confirmed by automated measure- 
ments. The properties of the selected Clumpy galaxies 
are consistent with the expectations for gravitationally 
unstable disks rather than mergers, in particular specific 
star formation rates (sSFR) on the Main Sequence of star 
formation rather than in a starbursting modtQ Velocity 
fields, available for some of these objects in the Litera- 
ture, are consistent with rotating disks. Thus, we have 
compared a sample of violently unstable disk galaxies, 
similar to many z > 1 star-forming galaxies, to a sample 
of more stable, spiral-like disk galaxies, with lower gas 
fraction and weaker instabilities. 

Utilizing the Mass-Excitation (MEx) diagnostic intro- 
duced by Juneau et al. (2011), we find that violently un- 
stable (Clumpy) systems frequently host an AGN, which 
is rarely the case in more stable galaxies. Comparisons 
of results from the MEx and from other narrow-line di- 
agnostics indicate that the much higher [Oiii]/H/3 ratios 
observed in clumpy systems cannot simply result from 
lower gas-phase metallicity, and actually require frequent 
Scyfert 2-type activity. X-ray signatures are in general 
too weak to be securely identified in individual systems 
- although one system has a hard X-ray detection and 
a few cases have a soft X-ray excess w.r.t star forma- 



tion. Stacking reveals a clear X-ray excess in clumpy disk 
galaxies, which significantly exceed the X-ray contribu- 
tion from star formation. This confirms the identification 
of frequent AGN in these systems. Intrinsic AGN bolo- 
metric luminosities are conservatively estimated to be of 
a few 10 43 ergs" 1 and potentially of a few 10 44 ergs" 1 , in 
Clumpy disks with stellar masses of a few 10 10 M Q at 
z^0.7. We also find that these AGN are, on aver- 
age, substantially obscured in X-rays, but not necessarily 
Compton thick. 

Clumpy disk galaxies appear to have higher gas frac- 
tions than stable disks, but by a factor lower than two in 
our mass-matched and redshift-matched samples. Thus, 
Clumpy disks contain somewhat larger gas reservoirs^, 
but not to the point of directly explaining the strong 
enhancement of the [Om] and X-ray luminosities by a 
factor of about five. An associated dynamical process 
for AGN feeding, triggered by the presence of clumps, 
is thus required. It may be either black hole growth in 
the giant clumps (Elmegrecn et al. 2008) and/or gas in- 
flows toward a central BH triggered by the presence of 
giant clumps - which can sometimes consist in clump mi- 
gration and central coalescence, but should be triggered 
even when the clumps are not merging centrally (Bll). 
We cannot specify the location of the AGN in our sample 
as we did not use spatially-resolved tracers. The average 
AGN luminosity is about five times higher in our Clumpy 
disk sample than in the Stable disk control sample, al- 
though this cannot directly constrain the AGN fraction 
among clumpy galaxies in general, since our sample was 
emission line-selected. Given the frequency of Clumpy 
disk instabilities at high redshift, and the typical AGN 
luminosity of at least few 10 43 ergs -1 in many of these 
systems, a large part of today's BH mass in galaxies with 
stellar masses of 10 10 " 11 M could be fueled in this pro- 
cess. 

Violent disk instabilities should persist longer in lower- 
mass galaxies before being replaced by weaker secular 
modes (b ars and spirals). This ha s been tentatively ob- 
served by lElmegreen et aTl ()2009al ). is supported by our 
own sample (see Section 3.4), and is also f ound in sam- 
ples of zoom-in cos molog ical simulations (|Martig et al.l 
I2012t iKraliic et all 120121) . This is naturally explained 
by th e evolution of cold acc retion with redshift and 
mass (|Dekel fc Birnboimll2006f) and of gas consumptio n 
with mass and metallicity (jKrumholz fc Dekell 120 111) . 
both of which event ually result in disk stabilization 
iCacciato et al.l (|2012f) . If violent disk instabilities in 
galaxies fed by cold gas accretion drive the bulk of the 
cosmic SF history, a downsizing in the termination of star 
formation should result, i.e. lower-mass galaxies con- 
tinue toforjnstarsa£tively at lower redshifts, as observed 
(e.g., Juneau et al.ll2005D. The downsizing of BH growth 
(e.g., 



Hasinger et al.ll2005t lLabita et aL1l26o9f) could also 



naturally result from the fueling by violent disk instabil- 
ities at high redshift in high-mass galaxies, and down to 
intermediate-redshift in lower-mass galaxies, as a natural 
consequence of the mass and redshift dependence of disk 
instabilities. 



alt hough a large fraction o f mergers are also on the Main Se- 
quence (Kartaltcpc et al. 20fJ). 

9 as expected to trigger their instability 
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APPENDIX 

We describe in the appendices the detailed methods used for stellar mass and emission line ratio measurements. 

STELLAR MASSES 

The stellar masses used in this work were derived from if -band luminosit ies (iRetzlaff et al1l2010f) . assuming a 
redshift-dependent mass-to-luminosity ratio for star-forming galaxies given by lArnouts et al.l ( 120071 equation 2). A 
corr ection factor of .2 dex was applied to convert stellar masses from a Salpeter IMF (as used in lArnouts "et~aT1l200l 
to a lChabrierl (|2003l ) IMF (as appropriate for the MEx diagnostic, Jll). 

We performed SED fitt ing for objects wi th reliable luminosities fro m the near-ultraviol e t to m id-infrared in the 
FIREWORKS catalogue (IWuvts et all 120081 This was done using the IBruzual fc Charlotl (|2003D stellar population 
models, and assuming exponentially-decaying single-burst star formation histories (so-called tau-models). We used 
stellar ages from 10 Myr to 13.5 Gyr, and a Calzetti extinction law with E(B-V) from 0.0 to 0.7. The SED fitting 
results yielded stellar masses lower by 0.2-0.3 dex than those obtained using a redshift-dependent M /Lk ratio (for 
both smooth and clumpy disks). This is in fact consistent with the findings by iMaraston et all ([2010!) on the effects 
of tau-models on stellar mass estimates for star-forming galaxies: using such star-formation histories is expected to 
undcr-estimat e the real stellar ma ss by 0.2-0.3 dex. For this reason, we preferre d to use a redshift-dependent M/Lx 
ratio. Indeed. fBitsakis et al.l(|2011| ) have shown, using the lda Cunha et all (|2008| ) SED models, that redshift-dependent 
M/L ratios is in good agreement with the results from SED fitting when the assumed star formation histories are not 
limited to tau-models. 

We note that a more recent calibration of redshift-dependent M/Lk ratios was proposed more recently bv lllbert et al.l 
(|2010T l , but based only on tau-mo dels and hence likely to under-estimate the stellar mass of star- forming galaxies . 
We thus used the calibration from lArnouts et al.l (|2007[) , which is based on composite star- formation histories and, in 
the 0.5 < z < 1.0 range, is calibrated on galaxies in a stellar mass range very close to that of our samples. 

We attribute a uniform uncertainty of 0.3 dex to all stellar mass estimates, because of the uncertainties on the M/Lk 
calibration (0.21 dex) and if-band luminosities (<0.15 dex). Note that even SED fitting with tau-m odels would yield 
result s within this uncertainty. Stellar masses were estimated independently for objects C7 and S5 bv lRodrigues et al.l 
(|2008ft . and the independent estimated from these authors are very similar to ours (within 0.1 dex, i.e., well within 
the error bar). 

We also checked consistency with dynamical masses in two cases, using velocity fields published in lYang et al.l (|2008| ). 
Object S3 has a rotation velocity V ~ 160kms _1 for an optical radius R ~ 14 kpc and optical inclination i ~ 45°, 
indicating the dynamical mas s M r \ vn = RV 2 /G within the optical radius. Assuming that two thirds of the mass within 
this radius is baryonic (e.g., iForster Schreiber et al.l 12006ft and that one third of the baryons in z^0.7 star-forming 
galaxies are in the form of gas ( Tacconi et al.l I2010T ) . this gives a stellar mass log(M */ M Q ) = 10.8. For object, C7, 

R ~ 8 kpc, V ~ 105km s _1 , i ~ 50°, which results in log(M */M ) = 10.1. These dynamical estimates (log(M*)=10.8 
and 10.1 for S3 and C7, respectively) are consistent with our measurements (log(M*)=10.9 and 10.2, respectively) 
within the assumed uncertainty of 0.3 dex. 

EMISSION LINE RATIOS 

We checked in the ESO database that the slit used during the observations includes at least 75% of the B-band 
emission for each selected system - in general, the whole galaxy is covered. We measured emission line fluxes by fitting 
a Gaussian function on lines. 

A linear continuum was estimated on a 80Arange on each side of the line, which provides a large enough range after 
excluding sky lines, other emission lines, and instrumental features. The flux was measured by integrating the fitted 
Gaussian on the wavelength range over which it exceeds 2.5cr, where a is the noise level estimated on the continuum. 
Given the measured noise level, we assumed a uniform uncertainty of 15% on measured fluxes and 0.2 dex on line 
ratios. Alternative flux measurement achieved by integrating lines over the range where the exceed 2.5cr without using 
a Gaussian fit give results that lie within this uncertainty. 

Underlying H/3 absorption could be robustly fitted in almost half of the cases. The resulting correction to the H/3 
flux was a factor ranging from 1.09 to 1.21, with a median value of 1.13. We applied a constant 1.20 to all spectra, 
which is conservative as it tends to (slightly) lower the AGN probabilities on the MEx diagnostic for Clumpy disks, as 
the average H/3 absorption is in fact stronger for Stable disks than for Clumpy disks (see Section l3~l~T) . US absorption 
can be measured on stacked spectra for clumpy and stable disks samples independently (see Section [3 .lj) : it implies an 
H/3 absorption correction factor of 1.16 for stable disks and 1.12 for clumpy disks, so we tend to slightly over-correct 
the H/3 flux for clumpy disks, but in a way that does not alter the results and anyway remains within the assumed 
uncertainty. These potential differences in H/3 correction between Clumpy and Stable disks, compared to the 1.20 
correction factor applied to all spectra in our analysis, correspond to an offset smaller than 0.05 dex on the MEx 
diagram, not changing significantly the estimated AGN probabi lities. 

We note that [Oiii]/H/3 and [Oil] /H/3 ratios were measured bv lRodrigues et all (|2008ft from a different spectroscopic 
survey for three of our targets, and are consistent with our own measurements and the assumed uncertainties. 

We acknowledge useful and insightful discussions with David Alexander, Damicn Le Borgne, Elisabete Da Cunha, 
Bruce Elmegrecn, Dcbra Elmegreen, Eric Emsellem, Natascha Forstcr-Schrcibcr, Dale Kocevski, Matthew Lehnert, 



19 



Mark Sarzi and Romain Teyssicr, and useful comments from an anonymous referee. We acknowledge support from 
the EC through grants ERC-StG-257720 (FB, SJ), the CosmoComp ITN (FB, SJ), grants ISF 6/08, GIF G-1052- 
104.7/2009, NSF AST-1010033 and a DIP grant (AD). 

REFERENCES 



Agertz, O., Teyssier, R., & Moore, B. 2009, MNRAS, 397, L64 

Alexander, D. M., et al. 2003, AJ, 126, 539 

Arnouts, S., et al. 2007, A&A, 476, 137 

Atlee, D. W., Martini, P., Assef, R. J., Kelson, D. D., & 

Mulchaey, J. S. 2011, ApJ, 729, 22 
Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5 
Balestra, I., et al. 2010, A&A, 512, A12 

Bargcr, A. J., Cowic, L. L., Bautz, M. W., et al. 2001, AJ, 122, 
2177 

Bitsakis, T., Charmandaris, V., da Cunha, E., et al. 2011, A&A, 
533, A142 

Boiler, T., Brandt, W. N., & Fink, H. 1996, A&A, 305, 53 
Bouche, N., et al. 2010, ApJ, 718, 1001 

Bournaud, F., Combes, F., & Semelin, B. 2005, MNRAS, 364, L18 
Bournaud, F., Dekel, A., Teyssier, R., Cacciato, M., Daddi, E., 

Juneau, S., & Shankar, F. 2011, ApJ, ApJ, 741, L33 (Bll) 
Bournaud, F., & Elmegreen, B. G. 2009, ApJ, 694, L158 
Bournaud, F., et al. 2008, A&A, 486, 741 
Brooks, A. M., Governato, F., Quinn, T., Brook, C. B., & 

Wadsley, J. 2009, ApJ, 694, 396 
Bruzual, G., & Chariot, S. 2003, MNRAS, 344, 1000 
Cacciato, M., Dekel, A., & Genel, S. 2012, MNRAS, 421, 818 
Ceverino, D., Dekel, A., & Bournaud, F. 2010, MNRAS, 404, 2151 
Chabrier, G. 2003, ApJ, 586, L133 
Chary, R., & Elbaz, D. 2001, ApJ, 556, 562 
Cisternas, M., et al. 2011, ApJ, 726, 57 

Comb es, F. 2001, Advanced Lectures on the Starburst-AGN , 

223, |aStro-ph/0010570l 
Cowie, L. L., Songaila, A., Hu, E. M., & Cohen, J. G. 1996, AJ, 

112, 839 

Cresci, G, Mannucci, F., Maiolino, R., et al. 2010, Nature, 467, 
811 

Cresci, G., Mannucci, F., Sommariva, V., et al. 2012, MNRAS, 
421, 262 

da Cunha, E., Chariot, S., & Elbaz, D. 2008, MNRAS, 388, 1595 
Daddi, E., Alexander, D. M., Dickinson, M., et al. 2007, ApJ, 
670, 173 

Daddi, E., Dannerbauer, H., Elbaz, D., Dickinson, M., Morrison, 

G., Stern, D., & Ravindranath, S. 2008, ApJ, 673, L21 
Daddi, E., Bournaud, F., Walter, F., et al. 2010, ApJ, 713, 686 
Daddi, E., Elbaz, D., Walter, F., et al. 2010b, ApJ, 714, L118 
Dekel, A., & Silk, J. 1986, ApJ, 303, 39 
Dekel, A., & Birnboim, Y. 2006, MNRAS, 368, 2 
Dekel, A., Sari, R., & Ceverino, D. 2009, ApJ, 703, 785 (DSC09) 
Dekel, A., et al. 2009, Nature, 457, 451 

Dickinson, M., & GOODS Team 2004, Bulletin of the American 

Astronomical Society, 36, 701 
Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 

604 

Di Matteo, P., Bournaud, F., Martig, M., Combes, F., Melchior, 

A.-L., & Semelin, B. 2008, A&A, 492, 31 
Di Matteo, T., et al. 2012, ApJ, 745, L29 

Dubois Y., Pichon, C, Haehnelt, M., et al. 2012, MNRAS, 3123, 

xxx l|arXiv:1112.2479|l 
Elbaz, D., & Cesarsky, C. J. 2003, Science, 300, 270 
Elbaz, D., Dickinson, M., Hwang, H. S., et al. 2011, A&A, 533, 

A119 

Ellison, S. L., Patton, D. R., Mendel, J. T., & Scudder, J. M. 

2011, MNRAS, 418, 2043 
Elmegreen, D. M., Elmegreen, B. G., & Hirst, A. C. 2004, ApJ, 

604, L21 

Elmegreen, B. G, & Elmegreen, D. M. 2005, ApJ, 627, 632 
Elmegreen, D. M., Elmegreen, B. G., Ravindranath, S., & Coe, 

D. A. 2007, ApJ, 658, 763 
Elmegreen, B. G, Bournaud, F., & Elmegreen, D. M. 2008, ApJ, 

684, 829 

Elmegreen, B. G, Bournaud, F., & Elmegreen, D. M. 2008b, 
ApJ, 688,67 

Elmegreen, B. G., Elmegreen, D. M., Fernandez, M. X., & 
Lemonias, J. J. 2009a, ApJ, 692, 12 



Elmegreen, D. M., Elmegreen, B. G., Marcus, M. T., Shahinyan, 
K., Yau, A., & Petersen, M. 2009b, ApJ, 701, 306 

Elmegreen, B. C, & Burkert, A. 2010, ApJ, 712, 294 

Emsellem, E., Greusard, D., Combes, F., Friedli, D., Leon, S., 
Pecontal, E., & Wozniak, H. 2001, A&A, 368, 52 

Fanidakis, N., Baugh, C. M., Benson, A. J., et al. 2011, MNRAS, 
410, 53 

Forster Schreiber, N. M., et al. 2006, ApJ, 645, 1062 
Forster Schreiber, N. M., et al. 2009, ApJ, 706, 1364 
Forster Schreiber, N. M., Shapley, A. E., Erb, D. K., Genzel, R., 

Steidel, C. C, Bouche, N., Cresci, G, & Davies, R. 2011, ApJ, 

731, 65 

Gabor, J. M., et al. 2009, ApJ, 691, 705 
Gammie, C. F. 2001, ApJ, 553, 174 

Genzel, R., Tacconi, L. J., Eisenhauer, F., et al. 2006, Nature, 
442, 786 

Genzel, R., et al. 2008, ApJ, 687, 59 
Genzel, R., et al. 2011, ApJ, 733, 101 

Georgakakis, A., Coil, A. L., Laird, E. S., et al. 2009, MNRAS, 
397, 623 

Giavalisco, M., et al. 2004, ApJ, 600, L93 

Grimm, H.-J., Gilfanov, M., & Sunyaev, R. 2003, MNRAS, 339, 
793 

Grogin, N. A., et al. 2005, ApJ, 627, L97 

Guo, Y., Giavalisco, M., Ferguson, H. C, Cassata, P., & 

Kockemocr, A. M. 2011, ApJ submitted. larXiv:1110 .3800 
Hasinger, G, Miyaji, T., & Schmidt, M. 2005, A&A, 441, 417 
Hcckman, T. M., Kauffmann, G., Brinchmann, J., Chariot, S., 

Tremonti, C, & White, S. D. M. 2004, ApJ, 613, 109 
Heckman, T. M., Ptak, A., Hornschemeier, A., & Kauffmann, G. 

2005, ApJ, 634, 161 
Hopkins, P. F., & Hernquist, L. 2006, ApJS, 166, 1 
Hopkins, P. F., Hernquist, L., Cox, T. J., Di Matteo, T., 

Robertson, B., & Springel, V. 2006, ApJS, 163, 1 
Ilbert, O., et al. 2010, ApJ, 709, 644 
Inoue, S., & Saitoh, T. R. 2012, MNRAS, 422, 1902 
Johansson, P. H., Burkert, A., & Naab, T. 2009, ApJ, 707, L184 
Juneau, S., et al. 2005, ApJ, 619, L135 

Juneau, S., Dickinson, M., Alexander, D. M., & Salim, S. 2011, 

ApJ, 736, 104 (Jll) 
Kauffmann, G. et al. 2003, MNRAS, 346, 1055 
Kannappan, S. J. 2004, ApJ, 611, L89 

Kartaltepe, J. S., Dickinson, M., A lexander, D. M., et al. 2012, 

ApJ submitted, larXiv:1110.4057l 
Kassin et al., ApJ submitted 

Keres, D., Katz, N., Weinberg, D. H., & Dave, R. 2005, MNRAS, 
363, 2 

Kocevski, D. D., Faber, S. M., Mozena, M., et al. 2012, ApJ, 744, 
148 

Kocevski, D., et al., in preparation 

Kraljic, K., Bournaud, F., & Martig, M., 2012, ApJ submitted 
Krumholz, M., & Burkert, A. 2010, ApJ, 724, 895 
Krumholz, M. R., & Dekel, A. 2010, MNRAS, 406, 112 
Kr umholz, M. R., fc Dekel, A. 2011, ApJ submitted 

HarXiv:1106.0301l l 
Labita, M., Decarli, R., Treves, A., & Falomo, R. 2009, MNRAS, 

399, 2099 
Lamareille, F. 2010, A&A, 509, A53 

Lehnert, M. D., Nesvadba, N. P. H., Tiran, L. L., Di Matteo, P., 
van Driel, W., Douglas, L. S., Chemin, L., & Bournaud, F. 
2009, ApJ, 699, 1660 

L 'Huillier, B., Com bes, F., & Semelin, S., A&A submitted, 
larXiv: 1108.42471 

McGaugh, S. S. 1991, ApJ, 380, 140 

Magorrian, J., et al. 1998, AJ, 115, 2285 

Maraston, C, Pforr, J., Renzini, A., Daddi, E., Dickinson, M., 
Cimatti, A., & Tonini, C. 2010, MNRAS, 407, 830 

Marconi, A., Risaliti, G, Gilli, R., Hunt, L. K., Maiolino, R., & 
Salvati, M. 2004, MNRAS, 351, 169 



20 



Martig, M., Bournaud, F., Teyssier, Ft., & Dekcl, A. 2009, ApJ, 
707, 250 

Martig, M., Bournaud, F., Croton, D. J., De kel, A., & Teyssier, 

R. 2012, ApJ submitted. larXiv:1201. 10791 
Martizzi, D., Teyssier, R., & Moore, B. 2012, MNRAS, 420, 2859 
Merritt, D., & Ferrarese, L. 2001, MNRAS, 320, L30 
Mortlock, D. J., Warren, S. J., Venemans, B. P., et al. 2011, 

Nature, 474, 616 
Mullancy, J. R., Alexander, D. M., Huynh, M., Goulding, A. D., 

& Frayer, D. 2010, MNRAS, 401, 995 
Mullaney, J. R., Pannella, M., Daddi, E., et al. 2012a, MNRAS, 

419, 95 

Mullaney, J. R., Daddi, E., Be thermin, M., et al. 2012b, ApJL 

submitted, larXiv:1204.2824l 
Neichel, B., et al. 2008, A&A, 484, 159 

Netzer, H., Mainieri, V., Rosati, P., & Trakhtenbrot, B. 2006, 

A&A, 453, 525 
Netzer, H. 2009, MNRAS, 399, 1907 
Noguchi, M. 1999, ApJ, 514, 77 

Nordon, R., Lutz, D., Genzel, R., et al. 2012, ApJ, 745, 182 
Popcsso, P., et al. 2009, A&A, 494, 443 
Puech, M. 2010, MNRAS, 406, 535 

Ranalli, P., Comastri, A., & Setti, G. 2003, A&A, 399, 39 
Ranalli, P., Comastri, A., Zamorani, G., et al. 2012, A&A in 

press, larXiv: 1204.44851 
Rctzlaff, J., Rosati, P., Dickinson, M., Vandame, B., Rite, C, 

Nonino, M., Cesarsky, C, & GOODS Team 2010, A&A, 511, 50 
Richstone, D., et al. 1998, Nature, 395, 14 
Riechers, D. A., Walter, F., Carilli, C. L., Bertoldi, F., & 

Momjian, E. 2008, ApJ, 686, L9 
Robertson, B. E., & Bullock, J. S. 2008, ApJ, 685, L27 
Rodighiero, G., Daddi, E., Baronchelli, I., et al. 2011, ApJ, 739, 

L40 

Rodrigues, M., et al. 2008, A&A, 492, 371 
Salmi, F., Daddi, E., Elbaz, D. et al. 2012, ApJ in press, 
arXiv:0490394 



Savaglio, S., Glazebrook, K., Le Borgne, D., et al. 2005, ApJ, 635, 
260 

Schawinski, K., Treister, E., Urry, C. M., Cardamone, C. N., 

Simmons, B., & Yi, S. K. 2011, ApJ, 727, L31 
Shapiro, K. L., et al. 2008, ApJ, 682, 231 

Shaplcy, A. E., Steidel, C. C, Erb, D. K., ct al. 2005, ApJ, 626, 
698 

Shcn, Y., Richards, G. T., Strauss, M. A., et al. 2011, ApJS, 194, 
45 

Shuder, J. M., & Osterbrock, D. E. 1981, ApJ, 250, 55 

Soltan, A. 1982, MNRAS, 200, 115 

Springel, V., & Hernquist, L. 2005, ApJ, 622, L9 

Tacconi, L. J., Genzel, R., Smail, I., ct al. 2008, ApJ, 680, 246 

Tacconi, L. J., Genzel, R., Neri, R., et al. 2010, Nature, 463, 781 

Teplitz, H. I., et al. 2011, AJ, 141, 1 

Treister, E., Urry, C. M., Schawinski, K., Cardamone, C. N., & 

Sanders, D. B. 2010, ApJ, 722, L238 
Treister, E., Schawinski, K., Volonteri, M., Natarajan, P., & 

Gawiser, E. 2011, Nature, 474, 356 
van Starkenburg, L., van der Werf, P. P., Franx, M., ct al. 2008, 

A&A, 488, 99 

Trouille, L., Barger, A. J., & Tremonti, C. 2011, ApJ, 742, 46 
Vanzella, E., et al. 2008, A&A, 478, 83 

Walter, F., Riechers, D., Cox, P., et al. 2009, Nature, 457, 699 

Worsley, M. A., et al. 2005, MNRAS, 357, 1281 

Wright, S. A., Larkin, J. E., Law, D. R., Steidel, C. C, Shaplcy, 

A. E., & Erb, D. K. 2009, ApJ, 699, 421 
Wright, S. A., Larkin, J. E., Graham, J. R., & Ma, C.-P. 2010, 

ApJ, 711, 1291 

Wuyts, S., Labbe, I., Schreiber, N. M. F., Franx, M., Rudnick, G., 

Brammer, G. B., & van Dokkum, P. G. 2008, ApJ, 682, 985 
Xue, Y. Q., Luo, B., Brandt, W. N., et al. 2011, ApJS, 195, 10 
Yan, R., Ho, L. C, Newman, J. A., ct al. 2011, ApJ, 728, 38 
Yang, Y., et al. 2008, A&A, 477, 789 



